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ABSTRACT 

We present F435W (B), F606W (Broad V), and F814W (Broad /) coronagraphic images of the debris disk 
around j3 Pictoris obtained with the Hubble Space Telescope's Advanced Camera for Surveys. These images 
provide the most photometrically accurate and morphologically detailed views of the disk between 30 and 
300 AU from the star ever recorded in scattered light. We confirm that the previously reported warp in the inner 
disk is a distinct secondary disk inclined by ^ 5° from the main disk. The projected spine of the secondary 
disk coincides with the isophotal inflections, or "butterfly asymmetry," previously seen at large distances from 
the star. We also confirm that the opposing extensions of the main disk have different position angles, but we 
find that this "wing-tilt asymmetry" is centered on the star rather than offset from it as previously reported. The 
main disk's northeast extension is linear from 80 to 250 AU, but the southwest extension is distinctly bowed 
with an amplitude of ^ 1 AU over the same region. Both extensions of the secondary disk appear linear, but 
not collinear, from 80 to 150 AU. Within ^ 120 AU of the star, the main disk is ^ 50% thinner than previously 
reported. The surface-brightness profiles along the spine of the main disk are fitted with four distinct radial 
power laws between 40 and 250 AU, while those of the secondary disk between 80 and 150 AU are fitted with 
single power laws. These discrepancies suggest that the two disks have different grain compositions or size 
distributions. The F606W/F435W and F814W/F435W flux ratios of the composite disk are nonuniform and 
asymmetric about both projected axes of the disk. The disk's northwest region appears 20-30% redder than 
its southeast region, which is inconsistent with the notion that forward scattering from the nearer northwest 
side of the disk should diminish with increasing wavelength. Within ^120 AU, the otf435w-'wf606w and 
'«F435W-w^F8i4w colors along the spine of the main disk are ^ 10% and ^ 20% redder, respectively, than those 
of (3 Pic. These colors increasingly redden beyond ^ 120 AU, becoming 25% and 40% redder, respectively, 
than the star at 250 AU. These measurements overrule previous determinations that the disk is composed of 
neutrally scattering grains. The change in color gradient at ^ 120 AU nearly coincides with the prominent 
inflection in the surface-brightness profile at 115 AU and the expected water-ice sublimation boundary. We 
compare the observed red colors within ^120 AU with the simulated colors of non-icy grains having a radial 
number density °^ and different compositions, porosities, and minimum grain sizes. The observed colors 
are consistent with those of compact or moderately porous grains of astronomical silicate and/or graphite with 
sizes > 0.15-0.20 pm, but the colors are inconsistent with the blue colors expected from grains with porosities 
> 90%. The increasingly red colors beyond the ice-sublimation zone may indicate the condensation of icy 
mantles on the refractory grains, or they may reflect an increasing minimum grain size caused by the cessation 
of cometary activity. 

Subject headings: circumstellar matter — planetary systems: formation — planetary systems: protoplanetary 
disks — stars: individual { j3 Pictoris) 
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1. INTRODUCTION 

Since the initial discoveries of cool (~ 100 K) du st around 
nearby stars by the Infrared Astronomical Satellite ( Aumann 
(1985) . (3 Pictoris has been the foremost example of a young, 
main-sequence star with a resolved circumstellar disk of dust. 
The disk likely comprises the debris from disintegrating bod- 
ies in a nascent planetary system rather than primordial dust 
from the dissipating protostellar nebula (Backman & Parescd 
[l993j .Artvmowicz 1997; Lagrange et al. 2000; Zuck ermanl 
120011) . Spectro scopic evidence of multitudinou s star- 
grazing comets ("Lagrange-Henr i et alJ il988; .Beust et alJ 
[1990; Vidal-Madiar et al. 1994, and references therein) has 
motivated models of the disk as an admixture of gas and 
dust from colliding and evaporating comets located within 
a few tens of AU from 3 Pic jLecavelier des Etan gs et al.l 
fT996t iBeust & Morbidellilll995 iThebault et al...200a . The 
cometary origin of the dust is supported by the detec- 
tion of broad, 10 /im silicate emission like that observed 
in the spectra of comets Halley, Kohoutek, and others 
Jjelesco & Knacke 1991; Knacke et al. 1993; Aitke n et alJ 
fl993>: .Weinberger et al..2003>:.Okamoto et al..2004 ). 

Ground-based, coronagraphic images of /3 Pic reveal 
an asymmetric, flared disk extending at least 1800 AU 
from the star and viewed nearly "edge-on" ( Smith & Terrile 
1^84; Paresce & Burrows 1987; Golimowski et al. 199^ 
Kalas & Jewitt.l995;.Mouillet et al..l997a:.Larwood & Kal^ 
2001h . High-resolution Hubble Space Telescope (HST) and 
adaptive-optics images show that the inner part of the disk 
20-100 AU from /? Pic) is warped in a manner consis- 
tent with the presence of a secondary disk that is inclined by 
~ 4° from the main disk and perhaps sustained by a massive 
planet in a similarly inclined, eccentric orbit (Burrow^_g^^ 
[1995';'Mou illet et alJ l997'F: Heap et al. 2000; Augere au et al. 
12001). HST and ground-based images also reveal concentra- 
tions of dust along the northeast extension of the disk about 
500-800 AU from the star that have been interpreted as an 
asymmetric system of rings formed, along with other asym- 
metries in the disk, after a close encounter with a passing star 
(jkalas et al. 2000, 2001; Larwood & Kalas 2001). Spatially 
resolved mid-infrared images show an asymmetri c inner disk 
having depl eted dust wit hin 40 AU of /3 Pic (Lagage & PantirJ 
[ip94; Panti n et alJl997|) and obhque clump s of emission 20- 
80 AU from the star ( W ahhai et alJl2003t IWeinberg er et ali 
12003; Telesco et al. 2005). These features suggest the pres- 
ence of noncoplanar dust rings whose locations conform to 
the mean-motion resonances of a putative planetary system. 

Constraints on the sizes of the dust grains observed 
in scattered light have been based upon multiband 
(BVRI) imaging studies of the disk in both unpolarized 
iParesce & Burrows 1987; Lecavelier des Etangs et al. 1993) 
and polarized ( Gledhill et al. 1991; Wolstencroft et al. 1995 ) 
light. The unpolarized images indicate that the disk is col- 
orless (within uncertainties of 20-30%) at distances 100- 
300 AU from (3 Pic, though its B-band brightness may be sup- 
pressed at 50 AU from the star.^^ This neutral scattering by 

96743 

Throughout this paper, we compute the projected dimensions of the 
disk in astronomical units (AU) using the trigonometric parallax of tt = 
0;'05187 ± 0"00051 (or I/tt = 19.28 ±0.19 pc) reported for (3 Pic by 
ICrifo et al. 1 1 997 ) based on astrometric measurements conducted with the 
Hipparcos satellite. Consequently, the projected distances reported in this pa- 
per may differ from those appearing in papers published before 1997, which 
were based on an erroneous distance of 16.4 pc to the star. 



the disk has been customarily viewed as evidence that the dust 
grains are much larger tha n the wavelengths o f the scattered 
light (> 1 ^m). However. IChini et alJ J 199 It) noted that the 
B-, V-, and /-band scattering efficiencies of silicate spheres 
were similar for grains with radii of 0.2-0.3 //m. Attempts 
to reconcile the neutral colors of the disk with the 10-25% 
polarizat ion of scattered light from the di sk ha ve been prob- 
lemati c. IVoshchinnikov & Kriigell Jl999h and iKrivova et alJ 
i2000h found that, although the polarization alone is best fit- 
ted with a grain size distribution with a lower limit of a few 
microns, the observed neutral colors can only be replicated 
by adding submicron-sized grains and lowering either the re- 
fractive index of the grains or the proportion of the smallest 
grains. Given the adjustments needed to match the polariza- 
tion models with the highly uncertain disk colors, a more pre- 
cise multicolor imaging study of the (3 Pic disk is warranted. 

In this paper, we present multiband coronagraphic images 
of (3 Pic's circumstellar disk obtained w ith HST's Advanced 
Came ra for Surveys (ACS) iFord et alj r2003; Gonzagaet aQ 
linos'). These images reveal the disk between 30 and 300 AU 
from the star with unprecedented spatial resolution, scattered- 
light suppression, and photometric precision. These qualities 
permit the measurement of the disk's optical colors with 3-10 
times better precision than previously reported from ground- 
based observations. By deconvolving the instrumental point- 
spread function (PSF) from each image, we accurately deter- 
mine the brightnesses, morphologies, and asymmetries of the 
two disk components associated with the warp in the inner 
disk. Our fully-processed images and results will likely serve 
as the empirical standards for subsequent scattered-light mod- 
els of the inner disk until the next generation of space-based 
coronagraphic imagers is deployed. 

2. OBSERVATIONS AND DATA PROCESSING 
2.1. ACS Imaging Strategy and Reduction 

Multiband, coronagraphic images of the A5v star [3 Pic 
were recorded on UT 2003 Oc tober 1 using the High Resolu- 
tion C hannel (HRC) of ACS (F ord et alJ2'003llGonzaga et alJ 
I2OO5I) . The HRC features a 1024 x 1024-pixel CCD detector 
whose pixels subtend an area of 0"028 x 0"025, providing 
a ^ 29" X 26" field of view (FOV). Beta Pic was acquired in 
the standard "peak-up" mode with the coronagraph assembly 
deployed in the focal plane of the aberrated beam. The star 
was then positioned behind the small (0"9 radius) occulting 
spot located approximately at the center of the FOV. HST was 
oriented so that the disk's midplane appeared approximately 
perpendicular to the 5" occulting finger and the large (1 "5 ra- 
dius) occulting spot that also lie in the FOV. Short, medium, 
and long exposures were recorded through the F435W (B), 
F606W (Broad V), and F814W (Broad /) filters over three 
consecutive HST orbits. All images were digitized using the 
default analog-to-digital conversion of 2 e~ DN~'. This se- 
quence of exposures was promptly repeated after rolling HST 
about the line of sight by ~ 10°. This offset changed the 
orientation of the disk in the FOV by ^ 10° and facilitated 
the discrimination of features associated with the disk from 
those intrinsic to the coronagraphic PSF. Immediately before 
the exposures of (3 Pic, coronagraphic images of the A7iv star 
a Pictoris were recorded through the same filters to provide 
reference images of a star having colors similar to those of 
(3 Pic but no known circumstellar dust. A log of all HRC ex- 
posures is given in Tabled 

The initial stages of image reduction (i.e., subtraction of 
bias and dark frames and division by a noncoronagraphic 
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flat field) were performed by the ACS image calibration 

S iipeline at the Space Telescope Science Institute (STScI) 
Pavlovskv et al.ll2b05l) . To correct the vignetting caused by 
the occulting spots, we divided the images by normalized 
"spot flats" that were appropriately registered to the approx- 
imate locations of the migrator y occu lting spots on the date 
of our observations (Krist et al.1 l2004l) . We then averaged the 
constituent images of each set of exposures listed in Table [2 
after interpolating over static bad pixels and eliminating tran- 
sient bad pixels with a conventional 3a rejection algorithm. 
We then normalized the averaged images to unit exposure 
time and replaced saturated pixels in the long-exposure im- 
ages with unsaturated pixels at corresponding locations in 
the shorter-exposure images. Throughout this process, we 
tracked the uncertainties associated with each image pixel. In 
this manner, we created cosmetically clean, high-contrast im- 
ages and meaningful error maps for each combination of star, 
filter, and roll angle. Figure[2shows 29" x 10" sections of the 
reduced F606W images of /3 Pic and a Pic obtained at each 
roll angle. 

2.2. Subtraction of the Coronagmphic PSF 

To distinguish the brightness and morphology of the disk 
from the diffracted and scattered light of j3 Pic, the occulted 
star's PSF must be removed from each image. By observing 
a Pic and (3 Pic in consecutive //ST orbits, we limited the dif- 
ferences between the coronagraphic PSFs of the two stars that 
would otherwise be caused by inconsistent redeployment of 
the coronagraph assembly, gradual migration of the occ ulting 
spot, or changes in HST's thermally driven focus cycles jRristl 
1^02). We measured the positions of the stars behind the oc- 
culting spot using the central peaks of the reduced corona- 
graphic PSFs (Figure^ that result from the reimaging of in- 
completely occulted, spherically aberrated starlight by ACS's 
corrective optics (iKrisBi200()i) . The positions of [3 Pic and 
a Pic differed by ^ 0.8 pixel 0"02). This offset causes 
differences between the coronagraphic PSFs that are large in 
the immediate vicinity of the occulting spot, but the residual 
light at larger field angles (> 5" from the star) after PSF sub- 
traction is ^ 10^'^ times fainter than the disk's midplane at 
those field angles (Krist 2000). 

Optimal subtraction of the coronagraphic PSF requires ac- 
curate normalization and registration of the filter images of 
the reference star a Pic with the corresponding images of 
/? Pic. Because direct images of the two stars were not ob- 
tained, we estimated the brightnesses of each star in each 
ACS bandpass using the HST synthetic photometry pack- 
age, Synphot, which h as bee n developed and distributed by 
STScI (Bushouse et a i] fT99l . In doing so, we used the op- 
tical spectra of the A5v stars 0' Serpentis and Praesepe 154 
(^unn & Strvker 1983) to approximate the spectrum of (i Pic. 
Likewise, we approximated the spectrum of a Pic with that 
of the A5iv star HD 165475B. These proxies yielded syn- 
thetic Johnson-Cousins photometry that closely match the 
Cousins BVRI measurements of a Pic and j3 Pic reported by 
[Bessel ( 1990). Assuming V magnitudes of 3.27 and 3.86 for 
a Pic and j3 Pic, respectively, we computed synthetic flux ra- 
tios, Fa/Fp, of 1.65, 1.75, and 1.90 for F435W, F606W, and 
F814W, respectively. We then divided the images of a Pic by 
these ratios to bring the integrated brightnesses of the refer- 
ence PSFs into conformity with those of (3 Pic. 

We aligned the normalized images of a Pic with the cor- 
responding images of j3 Pic using an interactive routine that 
permits orthogonal shifts of an image with subpixel resolu- 



tion and cubic convolution interpolation. The shift intervals 
and normalization factors (i.e., Fa/Fp) were progressively re- 
fined throughout the iterative process. We assessed the qual- 
ity of the normalization and registration by visually inspecting 
the difference image created after each shift or normalization 
adjustment. Convergence was reached when the subtraction 
residuals were visibly minimized and refinements of the shift 
interval or normalization factor had inconsequential effects. 
Based on these qualitative assessments, we estimate that the 
uncertainty of the registration along each axis is 0.125 pixel 
and the uncertainty of Fa/Fp in each bandpass is 2%. 

After subtracting the coronagraphic PSFs from each image, 
we transformed the images to correct the pronounced geo- 
metric distortion in the HRC image plane. In doing so, we 
used the coefficients of the biquartic-polynomial distortion 
map provided by STScI (Meurer et al. 2002) and cubic con- 
volution interpolation to conserve the imaged flux. We then 
combined the images obtained at each HST roll angle by ro- 
tating the images of the second group clockwise by 9°.7 (Ta- 
ble [Til, ali gning the respective pairs of images according to 
the previously measured stellar centroids, and averaging the 
image pairs after rejecting pixels that exceeded their local 3cr 
values. Again, we tracked the uncertainties associated with 
each stage of image processing to maintain a meaningful map 
of random pixel errors. 

We combined in quadrature the final random-error maps 
with estimates of the systematic errors caused by uncertainties 
in the normalization and registration of the reference PSFs. 
Other systematic errors from cyclic changes of HST's focus 
and differences between the field positions and broadband 
colors of a Pic and f3 Pic are negligible compared with the sur- 
face b rightness of [3 Pic's disk over most of the HRC's FOV 
JKrislll2000) . Our systematic-error maps represent the con- 
volved differences between the optimal PSF-subtracted image 
of (3 Pic and three nonoptimal ones generated by purpose- 
fully misaligning (along each axis) or misscaling the images 
of a Pic by amounts equal to our estimated uncertainties in 
PSF registration and Fa/Fp. The total systematic errors are 
1-5 times larger than the random errors within ^ 3" of (3 Pic, 
but they diminish to 10-25% of the random errors beyond 
^ 6" of the star We refer to the combined maps of random 
and systematic errors as total-error maps. 

Figure 121 shows the reduced and PSF-subtracted images of 
the disk in each ACS bandpass. Each image has been rotated 
so that the northeast extension of the disk is displayed hor- 
izontally to the left of each panel. The images have been 
divided by the brightness of /3 Pic in each bandpass derived 
from Synphot. The alternating light and dark bands near 
the occulting spot reflect imperfect PSF subtraction caused 
by the slightly mismatched colors and centroids of a Pic and 
(3 Pic. The bands perpendicular to the disk have amplitudes 
that are ^ 50-100% of the midplane surface brightnesses at 
similar distances from /3 Pic. These residuals preclude ac- 
curate photometry of the disk within 1"5 (^ 30 AU) of the 
star and anywhere along the direction of the occulting finger 
Along the midplane of the disk, the photometric uncertainties 
due to PSF subtraction are ^ 5-10% at a radius of r = 30 AU 
and less than 1 % for r > 60 AU. 

All calibrated surface brightnesses and colors presented in this paper are 
based upon the following Vega-based apparent magnitudes for /3 Pic obtained 
from Synphot: mF435w = 4.05, ffipeoew = 3.81, and mpsMW = 3.68. The sys- 
tematic zero-point errors are < 2% (Sirianni et al. 2005), and the estimated 
errors from imperfectly matched reference spectra are ~ 1-2%. 
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Figure|3lshows alternate views of the disk in which the ver- 
tical scale is expanded by a factor of four over that presented 
in Figure 12] and the vertical dimension of the disk's surface 
brightness is normalized by the brightness measured along 
the "spine" of the disk. (The spine comprises the vertical loca- 
tions of the maximum disk brightness measured along the hor- 
izontal axis of each image, after smoothing with a 3 x 3 pixel 
boxcar.) The expanded vertical scale exaggerates the warp 
in the inner disk first observed in images tak en with HST's 
Wide Field Planetary Camera 2 (WFPC2) by Burrows et ajj 
^995). The multiband images shown in Figures|2]and|3]may 
be directly compared with the unfiltered optical image of the 
disk obtained with the Space Telescope Imaging Spectrograph 
(S TIS) corona graph ( Gradv et al. 2003) and shown in Figure 8 
of lHeapetalJEbOOl) . 

2.3. Deconvolution of the "Off-Spot" PSF 

Accurate assessment of the chromatic dependencies of the 
disk's color and morphology requires the deconvolution of 
the unocculted instrumental PSF from each HRC filter im- 
age. This deconvolution of the "off-spot" PSF is especially 
important for the F814W images, because very red photons 
(A > 0.7 /im) passing through the HRC's CCD detector are 
scattered diffusely from the CCD substrate into a large halo 
that contributes significantly to the wide-angle component 
of the PSF ( Sirianni et al. 2005). Unfortunately, no collec- 
tion of empirical off-spot reference PSFs exists yet for the 
HRC coronagraph. Consequently, we can deconvolve the off- 
spot PSFs only approximately by using synthetic PSFs gener- 
ated by the Tiny Tim software package distributed by STScI 
dkrist & Hook 2004). Tiny Tim employs a simplistic model 
of the red halo that does not consider its known asymmetries 
dKrist et al.ll20()5bl) . but this model is sufficient for assessing 
the general impact of the red halo on our images of /3 Pic's 
disk. 

We generated model off-spot PSFs using the optical pre- 
scriptions, filter transmission curves, and sample A5v source 
spectrum incorporated in Tiny Tim. For simplicity, we ap- 
proximated the weakly field-dependent PSF in each band- 
pass with single model PSFs characteristic of the center of the 
FOV. The model PSFs extended to an angular radius of 10". 
We corrected the geometrically distorted model PSFs in the 
manner described in ^2.21 and then deconvolved them from 
the PSF-subtracted images of (3 Pic obtained at each roll an- 
gle. In doing so, we applied the Lucy-Richardson deconvolu- 
tion algorithm (Richardson 1972; Lucv 1974) to each image 
outside a circular region of radius 1 "5 centered on the sub- 
tracted star (The amplitudes and spatial frequencies of the 
PSF-subtraction residuals within this region were too large to 
yield credible deconvolved data.) The imaged FOV lacked 
any bright point-sources by which we could judge conver- 
gence of the deconvolution, so we terminated the computa- 
tion after 50 iterations. Examination of intermediate stages of 
the process showed no perceptible change in the deconvolved 
images after ^ 45 iterations. 

Figures |3 and |5] are the deconvolved counterparts of Fig- 
ures |2] and |3] respectively. The compensating effect of the 
deconvolution is especially evident when comparing the mor- 
phologies of the disk in Figures |3] and [S] The color-coded 
isophotes of the disk in each bandpass are much more similar 
in the deconvolved images than in the convolved images. The 
amplified, correlated noise in the deconvolved images is char- 
acteristic of the Lucy-Richardson algorithm when applied to 
faint, extended sources. It is a consequence of the algorithm's 



requirements of a low background signal, nonzero pixel val- 
ues, and flux conservation on both local and global scales. 
These requirements also account for the disappearance of the 
negative PSF-subtraction residuals near the occulting spot. 
The requirement of local flux conservation ensures reliable 
photometry in regions where the ratio of signal to noise (S/N) 
is large and where PSF-subtraction residuals are small, but it 
makes photometric measurements elsewhere less accurate and 
their uncertainties nonanalytic. 

3. IMAGE ANALYSIS 

The processed ACS/HRC images shown in Figures|2land|3] 
are the finest multiband, scattered-light images of /? Pic's in- 
ner disk obtained to date.^^ Earlier ground-based, scattered- 
light images show the usual effects of coarse spatial resolu- 
tion and PSF instability caused by variable atmospheric and 
local conditions ( Smith & Terrile 1984; Paresce & Burrov^ 
1987"; Golimowski et al. 1993; Lecavelier des Etangs et alj 
1993; Kalas & Jewitt 1995; Mouillet et al. 1997a b). the un- 
published BVRI WFPC2 images of the disk described by 
.Burrows et al. ( 1995) have comparatively low S/N ratios be- 
cause WFPC2 lacks a coronagraphic mode and directly im- 
aged starlight scatters irregularly along the surface of its CCD 
detectors. These conditions forced short exposure times to 
avoid excessive detector saturation and created irreproducible 
artifacts in the PSF-subtracted images. Moreover, the con- 
struction of a WFPC2 reference PSF from images of p Pic 
obtained at several roll angles allowed possible contamination 
of the reference PSF by the innermost re gion of the dis k. The 
unfiltered STIS images of the disk (Heapetal. 2000) com- 
pare favorably with our HRC images, notwithstanding their 
lack of chro matic information and partial pupil apodization 
JGradv et'aiji200 3). Both sets of images show the same region 
of the disk, though STIS's narrow occulting wedge permitted 
imaging of the disk's midplane about 0"4 (8 AU) closer to 
the star. The HRC's circular occulting spot and Lyot stop ex- 
posed the regions around the projected minor axis of the disk 
that were obscured in the STIS images by diffraction spikes 
and the shadow of the occulting bar. The HRC images have 
twice the spatial resolution of the STIS images, and they ex- 
hibit better S/N ratios in the regions of the disk between 150 
and 250 AU from the star 

3.1. Disk Morphology 

Figure |6l shows isophotal maps of our F606W images of 
the disk before and after deconvolution of the off-spot PSF. 
These maps are qualitatively similar to those generated from 
our F435W and F814W images. They confirm the morphol- 
ogy and asymmetries of the disk repo rted pr eviously by sev- 
eral grou£s_and described in detail bv iKalas & Jewitt ( 1995|) 
and iHeap et alJ ( 12000 ). The brightness asymmetry between 
the disk's opposing extensions is particularly evident along 
the spine beyond '^100 AU of f3 Pic. The asymmetric cur- 
vature of the isophotes on opposite sides of the spine and the 
inversion of this asymmetr y across the pr ojected minor axis of 
the disk are also apparent. iKalas & Jewitt (.1995) referred to 
this diametrical antisymmetry as the "butterfly asymmetry." 

Our images provide the first credible look at the sur- 
face brightness along the disk's projected minor axis. The 
convex isophotes in this region indicate that the pinched 
appearances of the disk in the images of .Smith & Terril^ 

Fully processed images and error maps in FITS format can be obtained 
via the World Wide Web at http://acs.pha.jhu.edu/~dag/betapic 
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(11984 ). lGoli mowski et al.l(ll9 93h.|LecaveIier des Etangs et alJ 
lll993l).lMo uilIet et al. (1997©, and, to a lesser degree, 
lKalas&.Tewitt.a995.) are artifacts of oversubtraction of the 
reference PSF and/or self-subtraction of the disk along its mi- 
nor axis. Other scattered-light images yielded no information 
in this region because of obscuration by co ronagranhic masks 
fraresce & Burrowi Il987t loiedhill et ^ I199U iHeap et all 
l200()i) . Our images are tinged by residuals from the im- 
perfect subtraction of the linear PSF feature seen in Fig- 
ure ^ but these residuals do not affect the overall contours 
of the isophotes. The isophotes are more widely spaced along 
the northwestern semiminor axis than along the southeastern 
semiminor axis. Such brightness asymmetry in optically thin 
disks is often attributed to the enhanced forward-scattering 
efficiency of the dust, which implies that the side of the disk 
closer to Earth is inclined slightly northwest from the line of 
sight to the star. This d eduction is consisten t with the in- 
clination determined by Kalas & JewittI ( Il995h from single- 
scattering models of the surface brightness along the projected 
major axis of the disk. 

3.1.1. The Main and Secondary Disks 

iMouillet et al.l(ll9971j) viewed the inner warp as a deformed 
and thickened region of the disk, perhaps caused and sus- 
tained by a planet in an inclined, eccentric orbit within 20 AU 
of f3 Pic. Heap et al. (2000) subsequently interpreted the warp 
as a blend of two separate disk components, the lesser of 
which is inclined from the main component by 4°.6. The lat- 
ter interpretation is supported by Figure which shows the 
ratios of our HRC images after and before deconvolution of 
the off-spot PSF. The division of the PSF-deconvolved filter 
images (Figure 0} by their corresponding convolved images 
(Figure|2} accentuates the small or narrow features in the disk 
that are most affected by blurring from the instrumental PSF. 
Thus the sharply-peaked midplane of the outer disk appears 
prominently in Figure The midplane of the inner warp is 
not as well defined as that of the outer disk, but the apparent 
separation of the midplanes beyond ^ 80 AU from the star 
indicates that the warp is a distinct secondary disk inclined 
from, and perhaps originating within, the main disk. 

Mimicking Heap et al. (2000), we determined the contribu- 
tions of the main and secondary disks to the composite verti- 
cal scattered-light profile (i.e., the scattered-light distribution 
perpendicular to the midplane of the disk) by fitting two sim- 
ilar, vertically symmetric profiles to the composite profile as 
a function of distance from the star. In doing so, we assigned 
a shape to the individual profiles by qualitatively assessing 
various analytic functions at a few locations along the disk 
in our PSF-deconvolved images. We found that a "hybrid- 
Lorentzian" function - two Lorentzians of different widths 
whose top and bottom parts, respectively, are smoothly joined 
at an arbitrary distance from their common centers - satisfac- 
torily matched each individual profile. (We ascribe no physi- 
cal significance to this hybrid-Lorentzian function, nor do we 
claim that it uniquely or optimally characterizes each profile.) 
We used a non-linear least-squares fitting algorithm based on 
the Levenberg-Marquardt technique (Press et al. 1992) to de- 
termine the positions, amplitudes, widths, and top-to-bottom 
transition zones of the main and secondary hybrid-Lorentzian 
profiles at distances of 30-250 AU from f3 Pic. The resulting 
fits of the composite vertical profiles were best in the regions 
80-130 AU from the star, where the two disk components are 
sufficiently bright and separated to allow unambiguous dis- 
crimination of both components. Figure |8] shows the best fits 



to the composite vertical profiles of both extensions of the 
disk observed in our F606W image at distances of 100 AU 
from the star 

Figure |9] shows the relative positions of the fitted compo- 
nents of the composite vertical profile as a function of distance 
from f3 Pic. For each filter image and disk extension, we per- 
formed linear least-squares fits to the traces of the main and 
secondary disks within the regions where they are accurately 
identified (80-250 AU for the main disk and approximately 
80-130 AU for the secondary disk). These fits are overplotted 
in Figure |9] The relative position angles of the main disk's 
northeast and southwest extensions (measured counterclock- 
wise from the horizontal axis) are about -0°.4 and 0°.5, re- 
spectively. The difference between these po sition angles is 
similar to the "wing-tilt asymmetry" noted bv lKalas & JewitJ 
( 1995) at distances of 100-450 AU from /3 Pic. However, our 
fitted lines appear to intersect at a point almost coincident with 
the star, rather than at a point in the southwest extension, as re- 
ported by Kalas & Jewitt. Expanding the vertical scale of our 
traces of the main disk in the F606W image (FigurefTol shows 
that the spine of the northeast extension conforms to the linear 
fit down to 80 AU, at which point the main and secondary 
disks can no longer be credibly distinguished. However, the 
spine of the southwest extension exhibits more curved than 
linear behavior over the range of the fit. This bow in the 
southwest extension appears in our F435W and F814W im- 
ages with equal clarity and therefore does not appear to be 
an artifact of our decomposition of the vertical scattered-light 
profile. 

Table|2llists the position angles of the linear fits along each 
extension of the secondary disk relative to the correspond- 
ing fits along the main disk. The average tilt of the sec- 
ondary disk's southwest extension matches the 4°. 6 derived 
by Hea p et alJ (jSoOO) for both extensions. However, we find 
that the average tilt of the northeast extension is ^ 25% larger 
than the values obtained by Heap et al. (2000) and by us for 
the southwest extension. Extrapolating the linear fits to the 
secondary disk outward (Figure|6|l shows that the convex in- 
flections of the butterfly asymmetry are aligned with the sec- 
ondary disk, which supports the notion that the asymmetry is 
caused by radiatively expelled dust from the secondary disk 
dAugereau e t al. 2001). On the other hand, extending the lin- 
ear fits to the secondary disk toward jS Pic (Figure O shows 
that its nonparallel extensions would not converge near the 
star if they extended sufficiently inward. Instead, the north- 
east extension would intersect the northeast extension of the 
main disk at a distance of ^ 30 AU from fi Pic. Corona- 
graphic images with a smaller occulting spot and higher pixel 
resolution are needed to determine if this intersection actually 
occurs. 

Although our analysis of the two- component struc ture of 
the disk qualitatively supports that of Heap et al. (2000), our 
results reveal some errors in the presentation and interpre- 
tation of their results. These errors stem from the incorrect 
portrayal of the secondary disk's tilt relative to the northeast 
and so uthwest ex tensions of the main disk in Figures 11- 
13 of Heap et al. (2000). In these figures, the plotted data 

- and consequently the relative slope of the secondary disk 

- are inverted about the midplane of the main disk. The ori- 
entations of the main and secondary disks should appear as 
the y do in Fig ures [THSlof this paper and in Figure s 8 and 10 
of Heap et al. (2000). Because of this inversion. Heap et alj 
concluded that the disk's southeast side was brighter than its 
northwest side and, consequently, that the southeast side was 
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inclined toward the line of sight. These conclusions are not 
supported by our isophotal m ap of the disk (Figu re |6|l or by 
the single-scattering models of lKalas & JewittI ( Il995i) . 

3. 1.2. Thickness of the Main Disk 

Previous scattered-light studies indicated that the com- 
posite disk is uniformly thick within 1 " (135 AU) of 
B Pic and flared beyond that region ( Smith & Terrile' 
1984; Paresce & Burrows 1987; Artvmowicz et al. 1990; 
GoHmowskiet al. 1993; Kalas & Jewitt 1995; Mouill et et al.1 
1997a: Heap et al. 2000). These characteristics are evident 
in Figures |3] and |5] Previous measurements of the disk's full 
width at half its midplane (FWHM) brightness vary signifi- 
cantly, however, which suggests that the disk's intrinsic pro- 
jected width is less than or comparable to the resolution of the 
images from which the measurements were made. Our de- 
convolution of the off-spot PSF from the HRC images - the 
first such effort exhibited for scattered-light images of (3 Pic's 
disk - permits us to assess this resolution dependency and to 
measure more accurately the intrinsic width of the disk as a 
function of d istance from (3 Pic. 

Figure nn shows the full widths of each extension of the 
composite disk at half and one tenth of the maximum bright- 
ness of the midplane at distances 20-250 AU from the star. 
These widths are shown in the panels labeled "FWHM" and 
"FWO.IM," respectively. The associated pairs of thin solid 
and dashed curves show the widths of the northeast and south- 
west extensions obtained from the F606W image before and 
after PSF deconvolution. The short, thick curves show the 
measured widths of the hybrid-Lorentzian curves that best 
fit the main component o f the c omposite disk in the region 
80-150 AU from /3 Pic (ELB- The FWHM and FWO.IM 
plots show progressive reductions in the measured width of 
the main disk after PSF deconvolution and decomposition of 
the main and secondary disks. The up-to-30% reduction in 
the FWHM after deconvolution of the narrow off-spot PSF 
(FWHM K, 0"05) shows that previously reported measure- 
ments, particularly those obtained from the ground, reflect 
more the angular resolution of the observations than the in- 
trinsic thickness of the disk. 

The deconvolved and decomposed curves show that the 
FWHM and FWO. IM of both disk extensions are nearly con- 
stant (~ 11-13 AU and 50-60 AU, respectively) within 
~ 120 AU of (3 Pic, and they increase approximately linearly 
with distance beyond 120 AU. This behavior is consistently 
seen in the F435W, F606W, and F814W images. The south- 
west extension appears to thicken more rapidly th an the north- 
east ex tension beyond ^160 AU, as noted bv .Kalas & JewittI 
Jl995h . but the correlated noise injected into the deconvolved 
images by the Lucy-Richardson algorithm (Figure |5} pre- 
cludes accurate quantification of the FWHMs from those im- 
ages. Reverting to the PSF-convolved images, we find that 
the FWHM and FWO.IM of the southwest extension are 
~ 10 AU and ~ 60 AU larger, respectively, than the corre- 
sponding values of northwest extension 250 AU from the star. 
Our finding is consis tent with the FWHM measurements of 
iKalas & JewittI J1995 ) when the different angular reso lutions 
of the two sets of observations are considered. Heap et al.l 
(|^00) reported no such thickness asymmetry from their STIS 
images. Their measurements of FWHM at discrete locations 
50-200 AU from j3 Pic are consistent with our PSF-convolved 
measurements of the narrower northeast extension. However, 
their FWO.IM measurements at these locations are 20-60% 
smaller than our PSF-convolved measurements at the same 



locations in each extension. These discrepancies may reflect 
differences between the S/N ratios of the HRC and STIS im- 
ages at scattered-light levels below ^ 10"^ of the stellar flux. 

3.1.3. Planetesimal Belts or Rings 

Reexamining earlier W FPC2 an d ground-based corona- 
graphic images of /3 Pic, IKalas et al. (2000) noted several 
localized density enhancements along the midplane of the 
disk's northeast extension at distances of 500-800 AU from 
the star They did not see this structure in the southwest ex- 
tension, so they interpreted the clumps as a system of eccen- 
tric, nested rings formed in the after math of a close enco unter 
between the disk and a passing star. IWahhai et alJ J2003 ) and 
Weinberger et al. ( 2003) later reported mid-infrared images of 
the disk that show a distinct warp within 20 AU of the star that 
is oppositely tilted from the secondary disk seen in Figure 
Wahhai et al. (2003) also noted clumps of emission in both 
disk extensions within 100 AU of /? Pic that are arranged in 
diametrically opposite pairs centered on the star They inter- 
preted these clumps as a series of noncoplanar rings caused 
by gravitational interactions of the disk with a planetary sys- 
tem. Telesco et al. ( 2005) did not observe the innermost warp 
in their mid-infrared images of the disk, but they attributed a 
bright clump of 12-18 //m emission at 52 AU in southwest ex- 
tension to a concentration of 0.1-0.2 /im silicate and organic 
refractory grains heated to ^ 190 K. High-spatial -resolution, 
mid-infrared spectra of the disk recorded at 3 AU intervals 
along the inner 30 AU of each extension also show con- 
centrated regions of e mission from 0.1-2 /im silicate grains 
("Okamoto et al."2004). The presence of such concentrations 
of small grains in the face of expulsive stellar-radiation pres- 
sure suggests that the grains are continuously created by col- 
lisions of planetesimals confined to rings or belts within the 
otherwise depleted inner region of the dust disk . 

The bright clumps detected by ' Kalas et all (|2£)00') in the 
northeast extension of the disk lie outside the FOV of 
our HRC images, and the innermost warp reported by 
IWahhai et all (12003') and' Weinberger et alJ ( 120031) is obscured 
by the HRC's occulting spot. However, the clumps of 
mid-infrared emission located 30-100 AU from the star 
(Wahhai etal. 2003; Telesco et al. 2005) are potentially ob- 
servable in our HRC images. Figure \^ is an isophotal 
map of the inner region of the disk obtained from our PSF- 
deconvolved F606 W image after smoothing with a 3 x 3 -pixel 
boxcar. The isophotal interval is 0.2 mag arcsec"^, which is 
similar to the interval between isophotes of 18 /im emission 
from the sa me region of the disk presented by Wahh ai et all 
( 12003.) and iTelesco et alJ 12005 ). The locations of the clumps 
of 18 u m emission ar e marked with the letters assigned to 
them bv IWahhai et al.' Our F606W isophotes show no evi- 
dence of such clumping in scattered optical light. Adjusting 
the isophotal interval and the smoothing factor does not alter 
this conclusion. Our images do not refute the existence of the 
clumps, however, because the scattered optical light and mid- 
infrared emission may emanate from different regions of the 
disk (£antin et al. 1997). 

3.2. Surface Brightness and Asymmetry 

Because f3 Pic's disk is viewed nearly "edge-on," its 
surface brightness has traditionally been parametrized with 
one or more power laws fitted along the midplanes 
(or spines) of its opposing extensions ( Smith & Terrild 
19841: lArtvmowicz et al 1989, 1990: Gledhill etal 199 iT 
Golimowski et alJ 119931: iLecavelier des Etangs et anfT99l 
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iKalas & JewitilQQStlMouillet et alJl997atlHeaD et al.l200(]l) . 

The number of fitted power laws and their indices have varied 
considerably, mostly because such fits are sensitive to errors 
in the stellar PSF subtraction. The innate stability of HST's in- 
strumental PSFs reduce these errors substantially, so we com- 
pare our quan titative measurements of the disk's spine (as de- 
fined in 42.21 with those obtained from the STIS images of 
iHeap et aL ti2000D . To do so, we assess the surface bright- 
nesses of the spines of both extensions from our ACS images 
both before and after PSF deconvolution. We initially ignore 
the two-component structure of the disk described in ^3.1.11 
because the two components cannot be distinguis hed by our 
profile-fitting algorithm within 80 AU of the star ( ^3.1.1> and 
because the main disk dominates the surface brightness be- 
yond 80 AU. We then compare our results for the composite 
disk with those obtained for its main and secondary compo- 
nents over the region in which they are credibly resolved. 

3.2.1. Horizontal Profiles of the Composite Disk 

Figure [Ol shows the logarithmic surface brightness pro- 
files measured along the spines of each extension of the 
composite disk before PSF deconvolution. The accompa- 
nying ±1(7 error profiles, which are derived from the total- 
error maps of the images (^^} and include ^ 3% uncer- 
tainty in the photometric calibration, show that these profiles 
are very accurate beyond ~ 1"5 30 AU) from the star 
The profiles clearly do not indicate a single power-law de- 
pendence of surface brightness with distance, as was deter- 
mined in so me of the earliest ground-based imaging studies 
of the disk (Smith & Terrile 1984; Artvmowicz et al. 19891 
[Gledhill et al. 1991; Lecavelier des Etanas et al. 1993). In- 
stead, the logarithmic profiles exhibit distinctly different lin- 
ear behavior at angular distances of 2"-3"5, 3"l-5"6, 6%- 
10", and 10"- 13 "4, with smooth transitions in slope be- 
tween these regions. Table |3] lists the slopes of the lin- 
ear least-squares fits to the logarithmic data (i.e., the in- 
dices of the radial power laws, r~", that best fit the surface- 
brightness profiles) within these four regions. The promi- 
nent changes in a at ~ 6" (^ 115 AU) from the star are 
well documented (Artvmowicz et al. 1990: Golimows ki et al.l 
[1993; Kalas & Jewitt 1995; Hear) et al. 2000), but the more 
subde changes at ~ 3"6 (~ 70 AU) and - 10" (- 195 AU) 
are noted here for the first time. Conversely, the pronounced 
decrease in the disk's ^-band surface brightness w i thin 7 "3 
(< 140 AU) reported bv' Lecavelier des Etangs et"al] lll993h is 
not observed in our F435W image. 

The prominent inflections in the profiles at ^ 115 AU are 
possible evidence for the sublimation of water ice from com- 
posite gr ains at dust temperatures greate r than 100-150 K 
jNakano" '1988'; 'Artvmowicz et al." 1990'; 'Golimowski et al.' 
U$93; Pantinetal. 1997; Li & Greenberg 1998). Applying 
their cometary dust model to (3 Pic's disk, iLi & Greenberg ! 
il99&) noted that the relative compositions of the grains 
should vary across several regions of the disk, as the inner 
and outer mantles of various ices and organic refractory ma- 
terial successively evaporate from, or condense upon, their 
silicate cores in accordance with the local dust temperature 
and stellar-radiation pressure. The more subtle changes in a 
at ^ 70 AU and ^195 AU may therefore mark the bound- 
aries over which such changes in composition (and conse- 
quently albedo) occur On the other hand, the decrease in 
a within ~ 70 AU may indicate a decrease in the number 
density of the gra i ns caused by an interior pl a netary system 
dLaeaee & PantinI 119941: iRoaues et all 119941: iLazzaro et all 



fT99l or, as we discuss in ^3.2.21 the unresolved superposi- 
tion of the main and secondary disks. 

The changes in a at angular distances of ^ 3 "6, 6", and 
10" from (3 Pic are also apparent in the midplane surface- 
brightness profiles extracted from the STIS images of the disk 
(Heapetal. 2000). To facilitate comparison with previous 
studies. Heap et al. fitted power-law functions to three regions 
of each extension. Two of these regions coincide approxi- 
mately with our regions 1 and 3 defined in Table |3j the third 
region bridges our regions 1 and 2. Our values of a in re- 
gion 1 (2"-3"5, or 39-67 AU, from (3 Pic) of both extensions 
differ by < 5% from those computed by Heap et al.l ll2000h for 
the positions of maximum flux along each extension. How- 
ever, our values of a in region 3 (6 "6- 10", or 12 7-193 AU, 
from f3 Pic) are ^ 10-15% smaller than those of Hea p et all 
for both extensions, i.e., our power-law fits in this region are 
significantly less steep than those determined from the STIS 
images. This discrepancy is puzzling, especially as this outer 
region of the disk is relatively insensitive to moderate errors 
in the PSF subtraction. 

To investigate the possible effect of the convolved instru- 
mental PSF on our power-law fits, we repeated this analy- 
sis on our PSF-deconvolved images in each bandpass. The 
surface-brightness profiles along the spines of each decon- 
volved disk extension are shown in Figure O The corre- 
sponding values of a for each region of the deconvolved im- 
ages are listed in Table |3] These values are significantly dif- 
ferent from their PSF-convolved counterparts only in region 4, 
where the S/N ratios are low and, consequently, the efficacy 
of the Lucy-Richardson algorithm is suspect. In region 3, 
the average difference between the associated convolved and 
deconvolved indices is ^ 1%, so it is unlikely that differ- 
ences between the HRC and STIS coronagraphic PSFs cause 
the ~ 10-15% d iscrepancies in the values of a reported by 
'Heap et al.' ("2000) and us for this region. 

Figure^ldisplays the PSF-deconvolved profiles in a man- 
ner that facilitates comparison of the two disk extensions in 
the F435W, F606W, and F814W images. The plots show the 
effects of the asymmetric values of a on the relative bright- 
nesses of the two extensions at equal distances from (3 Pic. 
The southwest extension is brighter than the northwest ex- 
tension in the region ~ 50-100 AU from the star, whereas 
the opposite condition exists for r > 150 AU. At r = 6" 
(~ 115 AU), the surface brightnesses of each extension are 
14.9, 14.6, and 14.4 mag arcsec'^ in F435W, F606W, and 
F814W, respectively. These values are 0.2-0.3 mag arcsec"^ 
brighter than the corresponding brightnesses measured from 
the PSF-convolved images. The latter measurements com- 
p are favorably with the gro und- based, j^-band measur ements 
of lGolimowski et alJ lll993h andlKalas & .Jewitt' ('19951) at r = 
6", but are ~ 1.5 mag arcsec"^ brighter than the and /- 
band measurements reported by Paresce & Burrows ( 19871) 
and Smith & Terrile ( .1984 ), respectively, for this fiducial dis- 
tance. The outer distance Umits of our photometric measure- 
ments (260 AU and 300 AU for the northeast and southwest 
extensions, respectively) are set by the HRC's FOV; they do 
not reflect an intrinsic asymmetry in the physical sizes of the 
extensions. 

Although the values of a within ^110 AU of f3 Pic (re- 
gions 1 and 2) are nearly constant across F435W, F606W, and 
F814W, those between 125 and 200 AU (region 3) decrease 
nonuniformly with increasing wavelength. These trends in- 
dicate that the optical colors of the disk are constant within 
^110 AU of the star, but they redden considerably between 



8 



GOLIMOWSKI ET AL. 



125 and 200 AU. The color gradient appears to flatten again 
from ^ 200 to 250 AU (region 4), but the relatively large un- 
certainties attached to a in this region preclude a definitive 
assessment. We confirm these trends in our subsequent anal- 
ysis of the disk's color ( ^3.3> . 

3.2.2. Horizontal Profiles of the Component Disks 

The resolution of the inner disk into two components 
prompts an assessment of the contributions of each compo- 
nent to the horizontal surface-brightness profiles of the com- 
posite disk. Figure[^shows the multiband, logarithmic pro- 
files of the main and secondary disks, extracted from the PSF- 
deconvolved images. The curves trace the maxima of the 
hybrid-Lorentzian profiles that best fit the vertical surface- 
brightness profiles of the two disks at 80-150 AU from (3 Pic 
( ^3. l.U . The secondary disk is ^ 1.5 mag arcsec"^ fainter 
than the main disk at 100 AU, as indicated in Figure |8l The 
profiles of the main disk show the same inflections at ^ 6" 
(r-^ 115 AU) exhibited by the composite disk, but the sec- 
ondary disk's profiles appear linear throughout this region. 
Whereas the relative brightnesses of the m ain disk's exten- 
sions invert at ^ 115 AU (as noted in ^3.2.1l for the composite 
disk), the southwest extension of the secondary disk is consis- 
tently brighter than its northeast extension from 80 to 150 AU. 

We fitted power laws to the profiles of both extensions of 
the main and secondary disks over the regions displayed in 
Figure^] taking care to avoid the inflection of the main disk 
at ^115 AU. Table |4] lists the indices, a, of these power 
laws. The main disk's indices at 80-108 AU are 15-35% 
smaller than those of the composite disk in the mostly over- 
lapping region 2 (defined in Table |3}, whereas its indices at 
127-150 AU are 5-20% larger than those of the composite 
disk in region 3. These indices are nearly constant across the 
B, V, and / bands at 80-108 AU, but they progressively de- 
crease with increasing wavelength at 127-150 AU. This be- 
havior mimics that of the composite disk in regions 2 and 3 
( ^3.2. 1> . which indicates that the main disk's colors progres- 
sively redden beyond the inflection at ^ 1 15 AU. Conversely, 
the secondary disk's indices progressively increase with in- 
creasing wavelength within the entire 80-150 AU region of 
the northeast extension, but they are effectively constant in 
the southwest extension. The former trend suggests that the 
colors in the northeast extension become bluer with increasing 
distance from (3 Pic. Altogether, these characteristics indicate 
that the grain populations and/or distributions in the main and 
secondary disks differ significantly at common projected dis- 
tances from the star. Moreover, if unresolved, the secondary 
disk can significantly affect models that use the scattered-light 
profiles to constrain the spatial distribution and composition 
of dust in the main disk. 

3.3. Optical Colors of the Disk 

Previous multiband-imaging studies of (3 Pic's disk showed 
that the optical colors of the star and disk are indistin- 
guishable within the studies' estimated photometric erro rs 
jParesce & Burrows 1987\ [Lecavelier des Eta ngs et al. 1993). 
These errors ranged from 20 to 30%, so actual color differ- 
ences of a few tenths of a magnitude would not have been 
deemed credible. The greater photometric precision of our 
ACS study enables a correspondingly more precise assess- 
ment of the disk's optical colors and scattering properties. We 
first examine the ACS colors of the composite disk to enable 
a comparison with the earlier ground-based studies. We then 



examine the colors of the main and secondary components de- 
rived from ou r deco mposition of the disk's vertical scattered- 
fight profile (EEB- 

3.3.1. Two-Band Flux Ratios of the Composite Disk 

Figure ini shows the F606W/F435W and F814W/F435W 
flux ratios of the disk obtained from the PSF-deconvolved im- 
ages shown in Figure 0] Because the images in Figure |3 are 
separately normalized by the stellar flux in each bandpass, the 
images in FigureFTTlrepresent the flux ratios of the disk rela- 
tive to those of the star Consequently, regions of the disk hav- 
ing a flux ratio of 1.0 exhibit the same otf435w-'Wf606w and 
'Wf435w-'«f814w colors (hereafter denoted F435W-F606W 
and F435W-F814W) as (3 Pic itself. Although the flux ra- 
tios along the projected minor axis of the disk are obscured 
by the residual artifacts of the PSF subtraction, those along 
the projected major axis are clear and uncorrupted. The im- 
ages show that the flux ratios are asymmetric about both pro- 
jected axes, with the largest ratios appearing in the north- 
northeast quadrant of the disk. Within 150 AU of the star, 
the F606W/F435W and F814W/F435W ratios on the south- 
east side of the disk are ^ 1.0-1.1 times those of the star, 
while the ratios on the northwest side are ~ 1.2-1.4 times 
those of the star The uncertainties of these ratios range from 
^ 5-10% along the spine to ~ 25% at projected vertical dis- 
tances of ±50 AU from the spine. Beyond 150 AU from 
the star, the color asymmetry between the two sides dimin- 
ishes substantially, but the uncertainties away from the spine 
are comparatively large (^ 25-75%). 

Fi gures l6l and 1171 show that the northwest side of the disk 
within 150 AU of (3 Pic is both brighter and redder than its 
southeast counterpart. If the brightness asymmetry indicates 
that the side of the disk nearer to Earth is tipped northwest- 
erly (§3.1), then the color asymmetry within ^ 150 AU indi- 
cates that the ensemble of dust grains in that region becomes 
increasingly forward-scattering with increasing wavelength. 
This behavior is inconsistent with the expected scattering 
properties of both compact and porous interste llar grains 
iDraine & Lee 1984; Wolff etal 1998: Voshchinnikov et alJ 
|2P05), so its cause is not evident. 

The vertical asymmetry of the composite disk's colors is 
strikingly well aligned with the spine of the disk, which ap- 
pears prominently in both panels of Figure As the S/N 
ratios of our flux measurements are largest along the spine, 
we hereafter restrict our analysis of the disk's colors to that 
region in order to assess their dependence on the horizontal 
projected distance from the star. Focusing on the spine also 
allows us to assess the impact of PSF convolution on the col- 
ors measured along this narrow feature of the composite disk. 

3.3.2. Colors along the Spine of the Composite Disk 

Figure [m shows the F435W-F606W and F435W-F814W 
colors of the composite disk relative to those of (3 Pic, before 
and after deconvolution of the off-spot PSF. We measured the 
colors along the spines of each disk extension after smooth- 
ing the F606W/F435 W and F8 14W/F435 W ratio images with 
a 7 X 7-pixel boxcar Although the PSF-convolved colors are 
less accurate than their deconvolved counterparts, their uncer- 
tainties are more suitably compared with those from previ- 
ous studies in which PSF deconvolution was not performed. 
The average uncertainties of both colors before PSF decon- 
volution are - 3% at 40-150 AU and ~ 8% at 150-250 AU, 
i.e., 3-10 times better than those obtained in previous studies 
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iParesce & Burrowsll987tlLecaveIier des Etangs et al.ll993l) . 
PSF deconvolution increases the differences between the 
F435W-F606W and F435W-F814W colors of the disk and 
star by ^ 0.03 and ~ 0.1 mag, respectively, at the cost 
of increased uncertainty. The relatively large increase in 
A(F435W-F814W) is due to the correction of blurring from 
both the instrumental PSF and the HRC's red-halo anomaly 
( ^2.3> . The general similarity of the respective solid and 
dashed curves in Figure shows that the color variations as 
a function of distance from [3 Pic are not significantly affected 
by PSF deconvolution. This condition is consistent with the 
practically invariant power-law indices that parametrize the 
midplane surface-brightness profiles in regions 1-3 of the 
composite disk before and after PSF deconvolution (Table|3}. 

The PSF-deconvolved curves in Figure show that the 
F435W-F606W and F435W-F814W colors of the disk are 
~ 0.1 mag and ^ 0.2 mag redder, respectively, than those 
of P Pic at distances of ^ 40-120 AU from the star. These 
color excesses are nearly constant within this region and have 
average uncertainties of ~ 8%. The F435W-F606W and 
F435W-F814W excesses respectively increase to ^ 0.2 mag 
and ^ 0.35 mag at 250 AU in the northeast extension and 
to ~ 0.25 mag and 0.4 mag at 250 AU in the southwest 
extension. The average uncertainties of both excesses be- 
yond 150 AU from the star ai-e - 23%, i.e., - 0.05-0.09 mag 
at 250 AU. The constancy of F435W-F606W and F435W- 
F814W within ^ 120 AU is consistent with the chromatically 
invariant power-law indices associated with regions 1 and 2 of 
the composite disk (Table O. Likewise, the steadily redden- 
ing colors from ^ 120-250 AU are consistent with the inverse 
proportionality of a and wavelength in region 3 of the disk. 

Our results contradict the longstanding notion that [3 Pic's 
disk scatters visible light neutrally and uniformly. Its red 
optical colors (relative to the star) are generally consis- 
tent with the colors of the disks surrounding the B9 Ve 
star HD 1415 69A (C lampin et al. 2003) and the G2 V star 
HD 107146 dArdila et al.i 2004. 2005), although those disks 
do not exhibit spatial color gradients like those seen in 
/3 Pic's disk beyond ^ 120 AU. The relatively red col- 
ors of j3 Pic's disk contrast markedly with the relatively 
blu e colo rs of the edge-on disk surrounding AU Microscopii 
jKrist etalJI2005al) . an M dwarf in a co-moving group of 
stars that includ es P ic (i Barrad o v Nava scues et aljll999t 
IZuckerman et an i2001l) . In ^4.2.61 we speculate that the dis- 
parate colors of the disks surrounding (3 Pic and AU Mic may 
be caused by small differences between the minimum sizes of 
dust grains in the two disks. 

3.3.3. Colors along the Spines of the Component Disks 

Figure [I9l shows the F435W-F606W and F435W-F814W 
colors of each extension of the main and secondary disks in 
the region 80-150 AU from j3 Pic. The colors are derived 
from the maxima of the hybrid-Lorentzian profiles that best 
fitted these components of the v ertical-scattered light profiles 
of the composite disk ( ^3.1.U . Because the colors depend 
upon the robustness of our fitting algorithm at each position 
along the spine of the composite disk, their uncertainties can- 
not be analytically determined. We estimate these uncertain- 
ties from the root-mean-squared (RMS) deviations computed 
for each disk extension. For the main disk, cr(F435w-F606W) ~ 
'^(F435w-F8i4w) ~ 0.02 over both extensions. For the sec- 
ondary disk, o'(F435w-F606w) ~ 0.1 over both extensions, and 
'''(F435W-F814W) ~ 0.1 and 0.2 over the southwest and north- 
east extensions, respectively. 



The colors of the main disk conform to those of the com- 
posite disk between 80 and 150 AU, as expected from its 
nearly fivefold superiority in brightness over the secondary 
disk. The colors of the secondary disk's southwest extension 
also match those of the main disk, but its northeast extension 
appears bluer than the main disk by ~ 0. 1-0.2 mag in F435W- 
F606W and ~ 0.1-0.4 mag in F435W-F814W. Moreover, the 
F435W-F814W color of this extension apparently becomes 
bluer with increasing distance from the star, as previously sug- 
g ested b y its steepening midplane surface-brightness gradient 
( ^3. 2. 21 and Table|3l. However, the relatively large uncertain- 
ties in this region of the disk preclude a definitive assessment 
of the colors and trends of the secondary disk. Within our esti- 
mated errors, the colors of the secondary disk are also consis- 
tent with being neutral and constant between 80 and 150 AU. 

4. DISCUSSION 

4.1. Impact of PSF deconvolution on observations and 
models 

Our analysis demonstrates that PSF deconvolution is an im- 
portant step toward accurately characterizing thin debris disks 
viewed edge-on or, more generally, less-inclined disks having 
narrow or compact structural elements. It is especially im- 
portant for multiband imaging studies designed to map the 
chromatic dependencies of a disk's morphology and surface 
brightness. For example, PSF deconvolution dispels the il- 
lusion that j3 Pic's disk thickens with increasing wavelength 
(Figures|3and|3} by restoring the displaced red flux to its orig- 
inal location along the disk's spine. This restoration of the true 
surface-brightness distribution contributes largely to our ob- 
servation that the main disk component is significantly redder 
than previously determined from ground-based measurement s 
iParesce & Burrows' 1987tlLecaveher des Etangs et alJl993ft . 
It also provides refined empirical constraints for past and fu- 
ture dynamical models that alternatively view the inner warp 
as a uniform, propagating deformation of a singl e disk or as 
two distinctly separat e, inclined disks ( Mouillet e t alJI 19973 
lAugereauet alJ2001l) . 

The improved measurements of the disk's morphology and 
surface brightness obtained from our PSF deconvolved im- 
ages provide a basis for improved scattered-light models of 
/? Pic's disk. However, deriving the spatial distribution and 
scattering properties of dust in an edge-on disk is compli- 
cated by the integrated effects of variations in both charac- 
teristics along the line of sight. This complication is com- 
pounded by the many asymmetries between the opposing 
extensions of (3 Pic's disk and by the existence of an in- 
clined secondary disk. Consequently, three-dimensional scat- 
tering models based on axisymmetric dust distributions - 
like those previo usly developed for the disks around (3 Pic 
j Artymowiczet aL 1989; .Kalas & Jewitt 1995) and AU Mic 
( Krist et al. 200li|) - have limited utility in the face of the 
morphological and photometric complexity of j3 Pic's disk re- 
vealed by our HRC images. 

Multicomponent, nonaxisymmetric models are clearly 
needed to reproduce the complicated scattered-light distribu- 
tion associated with [3 Pic's disk. Development of such mod- 
els is beyond the scope of this paper, but we identify several 
issues raised by our observations that should be addressed by 
future modelling efforts: 

1. The horizontal surface-brightness profiles of the main 
disk within ^150 AU are sensitive to contamination 
by flux from the secondary disk. Models of the dust 
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distribution within either disk should be based on de- 
convolved and deblended components of scattered light 
from both disks. 

2. The spines of the opposing extensions of the secondary 
disk are not collinear While the misaligned spines of 
the main disk (dubbed the "wing-tilt" asymmetry by 
iKalas & Jewitt 1995) are attributed to forward scatter- 
ing from a slightly inclined disk, those of the secondary 
disk are not as easily explained. The possible optical 
and/or dynamical causes of this phenomenon need to 
be explored. 

3. The dumpiness noted in mid-infrared images of 
the disk is not observed in optical scattered light, 
even though the grains believe d responsible for the 
clumps of mid-infrared emission (lOkamoto et alJ20'04l 
[Telesco et al. 2005) have sizes commensurate with ef- 
ficient scattering of optical light. This apparent incom- 
patibility may indicate that the scattered light and ther- 
mal emiss ion come from dif ferent regions along the line 
of sight ( Pantin et al.ll 19971) . If so, these disparate ob- 
servations may constrain the mean scattering properties 
of the grains. 

4. The optical colors of the composite disk are asymmet- 
ric about both projected axes of the disk. The gener- 
ally redder and brighter appearance of the northwest 
half of the disk is inconsistent with the expectation that 
the dust grains scatter red light more isotropically than 
blue light. This peculiar behavior may place interesting 
constraints on the sizes and compositions of the grains 
throughout the disk. 

5. The prominent changes in the horizontal surface- 
brightness profiles and color gradients at ^ 115 AU 
from P Pic suggest that these phenomena may have a 
common cause. The expected sublimation of water ice 
from micron-sized grains lyin g within ~ 100 AU of the 
star (iPantin et all 1 19971 lLi & Greenberg 1998) is one 
possible cause, but color variations caused by changes 
in the grain-size distribution must also be investigated. 

Although we defer these issues to future theoretical inves- 
tigations, we partly address the last issue here via a one- 
dimensional analysis of the colors observed along the spines 
of the component disks. 

4.2. Colors and Properties of the Dust 

iLi & Greenberd Jl998h successfully modeled the disk's 
continuum emission from near-infrared to millimeter wave- 
lengths and its 10 /im silicate emission feature by assuming 
that the dust is continually replenished by comets orbiting 
near, or falling onto, (3 Pic. In this scenario, the cometary 
dust grains are highly porous aggregates of primitive inter- 
stellar dust whose composition, molecular structure, and size 
and spatial distributions are altered by the stellar radiation 
enviro nments into which they are sputtered and transported. 
iLi & Greenberg (1998) did not apply this dust model to the 
disk's appearance in scattered light because the scattering 
properties of highly porous grains with appropriate compo- 
sitions had not been fully determ ined. Subsequent modeling 
of th ese properties ( Wolff etal. 1 998 : Vo shchinnikov et al.l 
l2005h now enables us to examine the porous-grain model 
from the perspective of the disk's optical colors. 



We begin by deriving a parameter that, in the absence of a 
detailed three-dimensional model of the disk, links the optical 
properties of a variety of grains with the colors observed along 
the spine of the disk. We apply this parameter first to the 
colors of the composite disk (Figure [Tst . as these colors are 
largely attributable to the main disk and their uncertainties 
are well established. Afterwards, we examine the less-certain 
colors measured for the secondary disk (Figure[T9l. 

4.2.1. Effective Scattering Efficiency of Midplane Grains 

The intensity of singly-scattered light measured along the 
midplane of an optically thin, edge-on disk is 

I{e)= j n{r)a,,^<^{e)Fa(^y^ dx, (1) 

where e is the angular distance from the star, n{r) is the num- 
ber density of grains at a distance r from the star, CTsca is the 
scattering cross section of grains of a given size and compo- 
sition, ^{9) is the phase function of the scattering angle 9, Fq 
is the stellar flux at an arbitrary radius ro, and x traverses the 
disk along the line of sight. The complicated and asymmetric 
surface brightness profiles of (3 Pic's disk suggest that n(r) 
cannot be expressed as a single analytic function throughout 
the disk. However, constraining n(r) within each region of 
each disk extension requires a complex modelling effort that 
is beyond the scope of this paper. 

For a less rigorous analysis of the disk's colors, we approxi- 
mate the surface brightness and number density profiles along 
the midplanes of both extensions with single radial power 
laws having fixed indices -a and -u, respectively. Because 
the ou ter radius of the disk is > 1800 AU (Larwood & Kal^ 
120011) . our lines of sight mostly traverse regions of the disk 
that lie well beyond the projected distances at which we see 
prominent inflections in the surface brightness profiles. The 
effective index a should therefore be biased toward the in- 
dices observed i n the outer regions of the disk, i.e., a w 4 
(Smith & Terril dl984 lKalas & Jewitt 1995). Nakano ( 199($ 
determined that v = a-1 if the observed scattered light orig- 
inates mostly from dust along those parts of line of sight 
near /3 Pic. This condition is valid unless the dust grains 
are strongly forward-scattering. Some models indicate that 
the sc attered light is dominated by such grains ( Pantin et al] 
Il997h . but the relationships between the sizes, compositions, 
and optical properties of the grains remain uncertain. Assum- 
ing for simplicity that the disk comprises grains that are at 
most moderately forward-scattering, we obtain 

«('")= • 

Adopting the disk geome try portrayed in Figure 1 of 
iBuitrago & Mediavillal(ll986l) . we rewrite Equation (1) as 

Fn nn ri fTT-eo 
^W = i;r^ / <Tsc. m^in'0d9, (3) 
sm e J 00 

where R is the distance to /3 Pic, and is flie scattering angle 
at the outer edge of the disk. The radius of the disk greatly 
exceeds the FOV of our images, so 9o « 0. The integral is 
then solely dependent upon the scattering properties of the 
grains, so we define the grains' effective scattering efficiency 
as 
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geff = / gsca H0) Sin'e dO, 

Jo 



(4) 



where Qsca = fsca/™^, and a is the grain size (radius). We 
use a common scattering phase function for small grains, 



m 



1- 



4-TT (l+^2_2^cos 6')3/2' 
where g is the scattering asymmetry parameter defined by 



(5) 



47r 



$(6*) cosOdn 



(6) 



and t/il is the unit solid angle (iHenvev & Greensteirilll94il) . 
Values of g — -I, 0, and 1 correspond respectively to back- 
ward, isotropic, and forward scattering. 

Substituting Equation (5) into Equation (4) and evaluating 
the integral, we obtain 



QeS 



(7) 



This result indicates that our view of the edge-on disk is pro- 
duced by mostly isotropic scatterers, which conforms to our 
initial assumption that the grains in the disk are not strongly 
forward-scattering. Its analytic form is a fortunate conse- 
quence of our choice of = 3, but the bias against forward- 
scattering grains is maintained for any choice of > 1 . Both 
2sca and g are functions of the dimensionless parameter x = 
2-Ka I A, where A is the wavelength of the scattered light. The 
mean Qeff for grains of size a, weighted by the stellar flux in 
a given bandpass, is 



(Geff) 



jFxTxdX ' 



(8) 



where F\ is the flux spectrum of /3 Pic and T\ is the filter 
transmission profile. We now compare (Qes) for grains of 
various compositions and porosities in the ACS B, Broad-V, 
and Broad-/ bands. 

4.2.2. Dust Grains within the Ice-Sublimation Zone 

Figure |20| shows (Geff) as a function of a for the F435 W, 
F606W, and F814W filters and five combinations of porosity 
and composition. None of these combinations includes icy 
mantles, so they are applicable only to the region of the disk 
within the ice sublimation limit (< 100 AU for micron-sized 
grains; Pantin et al. 1997; Li & Greenberg 1998). Figuresl^ 
and 120b depict (2eff) for compact (0% porosity) grains com- 
posed, respectively, of pure astronomical silicate ("astrosil") 
and equal parts astrosil and amorphous carbon (graphite). The 
curves are derived from tabulated values of Qsca and g for each 
composition based upon the work of Draine & Lee ( 1984). 
Figures l20t and I20fc show (Qeff) for 33% and 90% porous 
grains with equal amounts of astrosil and graphite inclusions, 
based on the values o f Qsca and g for such grains computed by 
IVoshchinnikov et aP (i2005.) . Figurel20tl depicts an intermedi- 
ate case of 60% porosity andj)ure astrosil derived from values 
of Qsca and g computed by Wolff etalJ (11998). For x > 25, 
we extrapolated the tabulated values of Qsca and g for porous 
grains to the geometric-optics limits of Qsca ^ 1 and g ^ 1 as 



era 

m 



As of this writing, tabulated values of Qsca and g for compact 
rains of astrosil and graphite are available on the World Wide Web at 
ltp://www.astro. princeton.edu/~draine/dust/dust. diel.html 



jc ^ oo. This extrapolation is acceptable because even the most 
porous grains are expected to be strongly forward-scattering 
forx ;2> 25 (N. Voshchinnikov, personal communication) and 
thus, for 1^ — 3, contribute negligibly to our edge-on view of 
the disk. 

Figures I20h and b represent the grain characteristics most 
often invoked when modelling scattered-light images of de- 
bris disks. When viewed edge-on, disks comprising com- 
pact grains with a ^ I exhibit similar scattering effi- 
ciencies in F435W, F606W, and F814W, i.e, they are neu- 
tral, albeit relatively inefficient, scatterers at optical wave- 
lengths. This condition is the basis of previous assertions 
that the reportedly neutral colors of /3 Pic's disk reflect a 
minimum grain size of several microns ( Pare see & Burrow^ 
11987; Lecavelier des Etangs et al. 1993). However, as noted 
by Chini et al. ( 1991), the values of (2eff) for smaller grains 
(a ss 0.2-0.3 iim) are also similar among the broad optical 
bands, so some submicron-sized grains are also neutral scat- 
terers at these wavelengths. Moreover, because small grains 
likely dominate the size distribution of grains in the disk 
(.Dohnanvi 1969; Li & Greenberg 1998), they more strongly 
influence the overall color of the disk than supermicron-sized 
grains. Nevertheless, the truly red colors of /3 Pic's composite 
disk (Figure I18> makes the debate over the cause of neutral 
colors in this case irrelevant. 

Figuresl20t. d, and e show that increasing porosity has two 
dramatic effects on (Qes). First, the pronounced peak at 0.2- 
0.4 /im diminishes rapidly for porosities > 60%. This phe- 
nomenon is caused by the similar behavior of Qsca for x < 10 
as porosity increases dWolff et al.lll998t IVosh chinniko v et all 
2005). Second, the peak broadens as it diminishes, caus- 
ing the ranges of predominantly blue-, neutral-, and red- 
scattering grains to broaden and shift to larger grain sizes. 
For 90% porosity, almost all grains with a <3 fim scatter 
Ught in F435W more efficiently than in F606W and F814W, 
and grains with a>3 /im are neutral scatterers. The values of 
(Geff) for blue-scattering grains with a < I /im are less than 
half those of the neutrally-scattering grains, but no distribu- 
tion of sizes for 90% porous grains will yield a disk with red 
colors. 

To compare the colors of an ensemble of grains having uni- 
form porosity and composition with the observed colors of 
f3 Pic's disk, we must determine for each bandpass the mean 
scattering cross section of the ensemble using (Geff) and a 
reasonable grain-size distribution. Applying this quantity to 
Equation (3), we express the intensity of scattered-light from 
the disk relative to the stellar flux in a given bandpass as. 



7Ta'{Qeff)^da, (9) 
da 



where = For^/R^ is the stellar flux measured at Earth, 
dno/da is the grain-size distribution, and amin and a^ax are 
the minimum and maximum sizes of the grains. The intrinsic 
color of the disk between bandpasses 1 and 2 is therefore. 



(10) 



OTi — m2 = — 2.5 log —r- 



Figure |^ shows the simulated F435W-F606W and 
F435W-F814W colors of the inner disk as functions of mini- 
mum grain size for the grain compositions and porosities pre- 
viously considered. To compute these colors, we assumed the 
values of (Qeff) from Figure|20|and a power-law size distribu- 
tion, dno °^ a~^-^ da, commonly attributed to dust produced 
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from planetesimal and particle collisions (lDohnanvilll969h . 
We computed the integral in Equation 9 for amm = 

0.01,0.02, 

5.00 /im and an arbitrary 

^max — 100 /im.^*^ Except for the 
case of 90% porosity, the observed colors along the spine of 
the composite disk within 120 AU (F435W-F606W « 0. 1 and 
F435W-F814W « 0.2; Figure [Tsl are simultaneously repro- 
duced by all combinations of porosity and composition when 
i^mm = 0.15-0.20 /im. The 90% porous grains yield only blue 
or neutral colors for any choice of Omin (as presaged by Fig- 
ure [20^), and therefore alone fail to explain the observed red 
colors of the inner disk. 

Our measured colors along the spine of the composite disk 
within 120 AU are, by themselves, insufficient to constrain 
tightly the composition and porosity of the grains in that re- 
gion. However, they do show that very porous grains do not 
contribute significantly to the integrated scattered light along 
the spine. This conclusion does not necessarily imply a lack 
of very porous grains, because the values of (Qeff) for such 
grains with a <2 /im are 10-20 times smaller than those of 
compact and moderately porous grains of similar size (Fig- 
ure|2()J. Thus, our results do not necessarily contradict those 
of Li & Greenberg (1998), who found that the dust models 
that best fitted the observed continuum and 10 fim silicate 
emission from the disk are those that feature extremely porous 
(> 95%) grains. It remains to be seen, however, whether a 
distribution of grain porosities can simultaneously satisfy the 
constraints provided by scattered-light and thermal images of 
the inner disk. 

4.2.3. Dust Grains beyond the Ice-Sublimation Zone 

Our estimate of flmin = 0.15-0.20 /im within 1 2 AU 
matches that obtained by Voshchinnikov & Kriig ell (Il999h 
using ground-based observations of the disk's colors and 
polarization beyond 115 AU ( ^aresce & Burrows 1987; 
iLecave lier des Etangs et al." "1993 ; Wolstenc roft et alJ 119951) 
and a dust model featuring an /?-band refractive index of niR = 
1.152-0.005/ and a grain-size distribution of dn °^ a~^-^ da. 
This refractive index is descriptive of both "dirty-ice" grains 
with 50% porosity and astrosil grains of 76% porosity. How- 
ever, our HRC images show that the disk steadily reddens 
beyond 120 AU, so the other dust models considered by 
[Voshchinnikov & Kriigel ( 1999) - and rejected because they 
produced colors that were too red compared with the nearly 
neutral colors observed from the ground - are in fact poten- 
tially viable models for this region of the disk. These mod- 
els feature larger refractive indices (corresponding to various 
combinations of compact and moderately porous grains of as- 
trosil, mixed or layered with dirty ice) and larger flmin- In an- 
other s tudy of the disk's polarization and colors, Krivova et al. 
(|2000) also favored models with moderate refractive indices, 
porosities < 50%, and depleted numbers of grains smaller 
than 2-3 /im. 

Although tabulated values of Qsca and g have not yet been 
published for moderately porous astrosil grains with icy man- 
tles, it is reasonable to expect that the optical colors of an 
edge-on disk of such grains vary with flmin in a manner sim- 
ilar to those depicted in Figures I21h-d. If so, then the in- 
creasingly red color gradient observed in j3 Pic's disk beyond 

We repeated this exercise for the porous-grain cases using the grain- 
mass distribution derived by Li & Greenberg 1 1998) from their thermal dust 
model that assumed 97.5% porosity and «(;-) « within the region 46 < 
r < 1 15 AU (Figure 1 Ic of Li & Greenberg 1998). The results were nearly 
identical to those obtained for the grain-size distribution o fiDohnanvi 119691) . 
so we discuss them no further. 



120 AU (Figure [Tst may reflect values of Omin that increase 
with distance from ~ 0. 15 um to perhaps 2 fim at 250 AU 
JVoshch innikov & Krugeri999^. A similar trend was pro- 
posed by Weinberger et al. (2003), based on the disappear- 
ance of the 10 /im silicate emission feature beyond 20 AU of 
the star. These phenomena are consistent with the scenario 
of a diminished presence of small (a < 2-3 /xm) grains near 
the ice-sublimation zone (> 100 AU), where cometary activ- 
ity ceases (Li & Greenberg 1998) and ra diation pressure set s 
«min to 1-10 /im, depending on porosity ( Artvmowicz 1988"). 
The steeper color gradients in the southwest extension of the 
disk suggest that flrnin increases more rapidly in that exten- 
sion, i.e., the number of small grains in the southwest exten- 
sion is smaller than in the northeast extension. This sugges- 
tion is consistent with the observed asymmetry between the 
polarizations of the two extensions ( Wolstencroft et al. 1995), 
which has been attributed to a 20-30 % larger pop ulation of 
small grains in the northeast extension jKrivova et al ..2000.) . 

4.2.4. Dust Grains in the Secondary Disk 

The largely uncertain colors of the secondary disk (' ^3.3.3> 
prevent us from rigorously comparing the grain characteristics 
of the main and secondary disks. However, the lack of inflec- 
tions in the secondary disk's color and surface-brightness pro- 
files (FigureO indicates that the grain populations of the two 
component disks are fundamentally different at 80-150 AU 
from the star If the inflections in the main disk's profiles 
at ~ 115 AU are solely caused by ice sublimation, then the 
composition and/or size distribution of the secondary disk's 
grains must be sufficiently different from those of the main 
disk that the ice-sublimation boundary in the secondary disk 
(if it exists) lies within or beyond the 80-150 AU range of 
our analysis. More precise photometry of the secondary disk 
is needed to determine whether the differences in the grains' 
composition or size distribution required for such a boundary 
sh ift are co mpatible with the colors of the disk. 

[Artvmowi cz ( 1997) argued that icy grains should not exist 
anywhere in (3 Pic's disk because they quickly photoevaporate 
in the face of the star's strong ultraviolet flux and because they 
are structurally brittle and unable to survive high-velocity col- 
lisions with other grains. If so, then other causes of the differ- 
ences between the color and surface-brightness profiles of the 
main and secondary disks must be considered. For example, 
the inflections in the main disk's profiles at ~ 115 AU may re- 
flect a sharp decrease in the number density of grains within 
that distance from the star. Lecavelier des Etangs et al. ( 199^ 
successfully modeled the inflection in the surface-brightness 
profile by imposing a lower limit on the sizes of comets or 
asteroids that travel from outer regions of the disk and re- 
tain enough volatile material to evaporate within 110 AU. If 
this scenario is correct, then the lack of inflections in the sec- 
ondary disk's profiles suggests that most of its dust is pro- 
duced by colliding or evaporating planetesimals orbiting near 
the star rather than by evaporating comets with large, ec- 
centric orbits. This interpretation is consistent with the hy- 
pothesis that the "warp" in the composite disk is caused by 
radiation-blown dust from colliding planetesimals that have 
been perturbed from the the innermos t part of the main disk 
by a giant planet in a n inclined orbit (iMouiUet etalJll997bt 
lAugereauetalJ2001h. 

4.2.5. Possible Effect of Irradiation on Disk Colors 

Recent photometric studies of Kuiper Belt Objects (KBOs) 
reveal that these objects exhibit a broad range of optical and 



MULTIBAND IMAGES OF DISK AROUND /3 PICTORIS 



13 



near- infrared colors (Lm i & Jewitlll l996t Ijewitt & Luul l 19981 
2(K)!llTegler & R omanish inlll998i l2003: Ba rucci et alJEOOlt 
Delsanti et alJl2004.) . Growing (but controversial) evidence 
suggests that KBOs are divided by their perihelia, q, into two 
color populations: KBOs with <7 < 40 AU have neutral-to- 
red intrinsic colors, and KBOs with ^ > 40 AU have only 
red intrinsic colors (0.5 < B-R 5, 0.75, relative to the Sun; 
iTegler & Romanishin|[T998l 120031 IPelsanti et al.ll200l . A 
possible cause of this color trend is the diminution of col- 
Usions and/or cometary activity in the outer Kuiper Belt 
dLuu & Jewitt 1996; Stern 2002; Dels anti et al.. .2004). This 
hypothesis rests on the notion that the "dirty-ice" surfaces of 
KBOs are progressively polymerized and reddened by high- 
energy solar and cosmic radiation unless they are recoated 
with primordial dust ejected by collisions between KBOs or 
the sublimation of H2O, CO2, and CO ices. 

The observed reddening of /3 Pic's disk beyond its ice- 
sublimation zone suggests a possible connection with the red- 
dening of distant KBOs. Could the reddening of f3 Pic's disk 
beyond ^115 AU be the result of irradiative polymerization 
rather than an increasing minimum grain size? Luu & JewitJ 
((1996) simulated the rate of irradiative reddening of KBOs 
using an inverse exponential function with an e-folding time 
of ~ 10^ yr Because the radiation is ov erwhelmingly solar 
within the heliopause (Gil-Hutton 2002'), we crudely apply 
this function to the region of (3 Pic's disk in our FOV by as- 
suming that the e-folding time scales with luminosity and the 
distance from the star In this situation, the B-V color of long- 
lived icy grains orbiting (3 Pic (L/Lq = 8.7; Crifo et al. 1997) 
would increase by ^0.3 mag in a few times 10^ yr. The life- 
times of grains not quickly expelled from the disk by radiation 
pressure are limited by either Poynting-Robertson drag or in- 
tergrain collisions, depending on their size and location within 
the disk. At 120-250 AU, the minimum lifetimes of grains 
larger than ~ 0.2 iim are constrained by collisions to ^ 0.2- 
2 Myr (Backman & Paresce 1993). Thus, the lifetimes of the 
grains observed in our HRC images are much too short for the 
observed color gradient to be caused by irradiative reddening. 
This condition likely persists throughout the outer disk, be- 
cause the grain lifetimes beyond ~ 500 AU (which are lim- 
ited by Poynting-Robertson drag) are ~ 1 % of the irradiative 
e-folding times at those distances. 

4.2.6. Comparing the disks around (3 Pic and AU Mic 

The disk around the M dwarf AU Mic, whose spa- 
tial velo city and age are similar to those of /? Pic 
jBarrado v Navascues et alJll999t) . is the only kn own debris 
disk exhibiting blue colors at optical wavelengths (iKrist et alJ 
|2005a). Moreover, the disk's F435W-F814W color becomes 
bluer with increasing distance from the star. The disk ex- 
tends at least 210 AU from the star and is viewed almost 
edge-on dKalas et aL 2004), so the predicted F435W-F814W 
colors shown in Figure 1211 can be compared to the observed 
ACS colors of the disk if the grains are not icy. I n this case, 
the measured F435W-F814W = -0.3 at 30 AU JKrist et al.1 
|2p05a) from AU Mic suggests that the disk may be com- 
posed of astrosil and/or graphite grains with fljnin < 0.1 fim 
for porosities < 60% or Oniin ~ 0.2 /im for 90% porosity. 
The possibility that AU Mic's disk contains smaller or more 
porous grains than are evident in f3 Pic's inner disk ( ^4.2.2> 
is consistent with the much lower radiation pressure exerted 
on suc h grains by AU Mic than by /? Pic (Kalas et al. 2004^ 
iKrist et al. 2005a ). On the other hand, the measured F435W- 
F814W = -0.5 at 60 AU from the star cannot be produced by 



any combination of grain composition and porosity shown in 
Figure|21 so the grains in AU Mic's disk may have a different 
mineral composition altogether 

5. SUMMARY AND CONCLUDING REMARKS 

We have presented B-, Broad V-, and Broad /-band corona- 
graphic images of the dusty debris disk around (3 Pictoris ob- 
tained with the High Resolution Channel of HST's Advanced 
Camera for Surveys. We have exploited the HRC's image 
resolution and stability by subtracting a well-matched coron- 
agraphic reference image and deconvolving the instrumental 
PSF from the multiband images. The resultant images pro- 
vide the most morphologically detailed and photometrically 
accurate views of the disk between 30 and 300 AU from the 
star obtained to date. 

Our PSF-deconvolved images confirm that the appar- 
ent warp in the disk < 100 AU from the star, which 
was previously observed by iBurrows et al\ (| 1995ft and 
Heapetal. (2003) a nd mo delled by Mouill et et alJ jl997bh 
and lAugereau et al. ( 2001), is a distinct secondary disk in- 
clined to the main disk by ^5°. The opposing extensions of 
the secondary disk are not collinear, but their outwardly pro- 
jected midplanes (spines) are coincident with the isophotal 
inflections previously seen at large distances and com monly 
called the "butterfly asymmetry" (Kalas & Jewitt 1995). The 
surface brightness profiles along the spines of the secondary 
disk from 80 to 150 AU can be fit with single power-law func- 
tions having indices of -3.7 < -a < -5.0. The lack of inflec- 
tions in the surface brightness profiles around the expected 
ice-sublimation boundary (^ 100 AU) suggests that the com- 
position and/or size distribution of grains in the secondary 
disk is different from those of the main disk. Altogether, these 
phenomena support the notion that the secondary disk and the 
butterfly asymmetry comprise radiatively expelled dust from 
colliding planetesimals in inclined orbits near the star rather 
than from evaporating comets in large, eccentric orbits. 

We confirm the "wing-tilt asymmetry" between the oppos- 
ing extensions of the main disk, but we find that the effect 
is centered on t he star rather th an offset toward the south- 
west extension (Kalas & Jewitt"1995). While the spine of 
the northeast extension appears linear 80-250 AU from the 
star, the southwest spine is distinctly bowed with an ampli- 
tude of ~ 1 AU. The vertical width of the main disk within 
~ 120 AU is nearly constant and is up to 50% narrower than 
previously reported. The clumpy structures observed in mid- 
infr ared images and spectra o£the disk ( Wahhai et al. 20(33 
[Weinberger et alJ 120031: lOkamoto et ali 2004; T elesco etall 
are not seen in our optical scattered-light images. The 
surface-brightness profiles along the spines of the main disk's 
extensions can be fit by four distinct power laws separated by 
inflections at 70, 117, and 193 AU. The power-law indices, 
a, match those of Heap et al. (2000) wel l within 70 AU of 
P Pic, but are 10-15% smafler than those of lHeap et al]tt2000) 
at 127-193 AU. This discrepancy cannot be attributed to dif- 
ferences in the instrumental PSFs. The indices within 150 AU 
change significantly after removing the contribution from the 
secondary disk, so the superposed effects of both disks must 
be considered in future models of (3 Pic's circumstellar dust 
distribution. 

The two-dimensional F606W/F435W and F814W/F435W 
flux ratios of the composite disk are nonuniform and asym- 
metric about the projected major and minor axes of the disk. 
Within 150 AU of $ Pic, the ratios on the southeast side of 
the disk are ~ 1.0-1.1 times those of the star, while those 
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on the northwest side are ^ 1.2-1.4 times those of the star. 
The redder appearance of the nearer northwest side of the 
disk is inconsistent with the expectation that forward scatter- 
ing from disk grains should diminish with increasing wave- 
length. The F435 W-F606W and F435W-F8 14W colors along 
the spine of the disk are ~ 0.1 mag and ^0.2 mag red- 
der, respectively, than those of (3 Pic at ^ 40-120 AU from 
the star. These color excesses increase steadily beyond ^ 
120 AU, respectively reaching ^ 0.2-0.35 mag and ^ 0.25- 
0.4 mag at 250 AU. These results contradict the longstand- 
ing notion that the disk consists of neutrally scattering grains 
with sizes larger than several iim JParesce & Burrowslll987t 
iLecavelier des Etan^s et al.ll993l) . 

We have compared the colors measured along the spine of 
the composite disk with those expected for non-icy grains 
having a number density and different compositions, 

porosities, and minimum grain sizes. We find that the ob- 
served F435W-F606W and F435W-F814W colors within the 
ice-sublimation zone (< 100 AU) are consistent with those of 
compact or moderately porous (P < 60%) grains of astronom- 
ical silicate and/or graphite with minimum sizes of ^ 0.15- 
0.20 fim. The observed colors are inconsistent with the blue 
colors expected from the very porous grains (P > 90%) that 
best reproduce the 10 fxm s ilicate emission featur e observed 
within ^ 35 AU of the star (■L i & Greenbergll998 '). The red- 
dening colors beyond ^120 AU may reflect the formation 
of "dirty ice" grains or an increasing minimum grain size be- 
yond the ice-sublimation boundary. The latter condition is 
consistent with the decreased production of submicron grains 
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as cometary activity diminishes. It is unlikely that the redden- 
ing of disk beyond ~ 120 AU is caused by uradiative poly- 
merization (as has been postulated for the reddest and most 
distant Kuiper Belt Objects) because the required irradiation 
time is hundreds of times longer than the expected lifetimes 
of the grains at those distances. 

Our ACS/HRC coronagraphic images of /3 Pic's disk are 
the finest multiband, scattered-light images of its inner region 
(30-300 AU) ever recorded. These images will not be su- 
perseded by the Terrestrial Planet Finder or other proposed 
extrasolar-planet imaging missions because the fields of view 
of those missions will lie well within the region obscured by 
the HRC's occulting spot. Comparable infrared images of the 
disk are expected from the James Webb Space Telescope if the 
current specifications for its coronagraphic-imaging modes 
are maintained. Thus, our observations and derived results 
should be standard references for comparative and theoreti- 
cal studies of circumstellar debris disks for at least the next 
decade. 
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TABLE 1 
Log of ACS/HRC exposures" 



Roll angle*" Exposures 



Star 


(deg) 


Orbit 


Filter 


No. 


X 


Time (s) 


ct Pic 


87.2 


1 


F435W 


1 


X 


40 










2 


X 


100 










2 


X 


225 








F606W 


1 


X 


15 










2 


X 


50 










2 


X 


200 








F814W 


1 


X 


40 










2 


X 


100 










2 


X 


200 


/3Pic 


110.0 


2 


F606W 


1 


X 


10 










2 


X 


100 










6 


X 


338 






3 


F435W 


1 


X 


10 










2 


X 


250 










3 


X 


653 






4 


F814W 


1 


X 


10 










2 


X 


250 










3 


X 


645 


/3Pic 


100.3 


5 


F606W 


1 


X 


10 










2 


X 


100 










6 


X 


338 






6 


F435W 


1 


X 


10 










2 


X 


250 










3 


X 


653 






7 


F814W 


1 


X 


10 










2 


X 


250 










3 


X 


645 



"Exposures recorded over seven consecutive HST orbits on UT 2003 October 1 . 
''Defined as the position angle of the y-axis of the raw HRC image, measured east 
of north. 



TABLE 2 

Relative position angles of secondary disk" 



Extension 


F435W 


F606W 


F814W 


Average 


Northeast 


6°.3 


5°.3 


6°. 3 


5°.9±0°.6 


Southwest 


-4°.4 


-4°.9 


-4°. 8 


-4°.7±0°.3 



"The position angles pertain to regions of the secondary disk between the vertical dotted lines 
in Figure|9] and are measured relative to the corresponding extensions of the main disk. 
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TABLE 3 

Power-law fits to surface-brightness profiles of composite disk 





PSF 






Power-law index (a) 




Extension 


Decon? 


Region 


F435W 


F606W 


F814W 


Northeast 


No 


1 


1.34 ±0.02 


1.34 ±0.02 


1.30 ±0.02 






2 


J 91 i 02 


1 Q 1 + 09 


1 QO + 09 






3 


4. 70_i_r) 07 


4 1 Q ± 09 


4 04 ± 09 






4 






J .^y ZZLL u.uo 




Yes 


1 


1.45 ±0.04 


1.43 ±0.03 


1.38±0.03 






2 


1.84 ±0.04 


1.83 ±0.03 


1.80 ±0.03 






3 


4.20 ±0.05 


4. 19 ±0.04 


4. 12 ±0.04 






4'-' 


2.80±0.15 


3.34±0.11 


3.34±0.12 


Southwest 


No 


1 


1.02 ±0.02 


1.06 ±0.02 


1.15±0.02 






2 


2.03 ±0.02 


2.03 ±0.02 


2.05 ±0.02 






3 


4.89 ±0.03 


4.76 ±0.02 


4.56 ±0.02 






4 


3.99 ±0.09 


4.00 ±0.05 


3.84 ±0.07 




Yes 


1 


1.08 ±0.04 


1.10 ±0.03 


1.21 ±0.04 






2 


1.93 ±0.04 


1.94 ±0.03 


1.97 ±0.03 






3 


4.82 ±0.06 


4.75 ±0.05 


4.60 ±0.05 






4-^ 


3.11±0.18 


3.53±0.12 


3.26±0.14 



'a is the index of the radial power-law, r^", that best fits the surface biightness of the spine within a given 
region of each disk extension. Equivalently, -a is the slope of the linear least-squares fit to the given region of the 
logarithm of the profile, as shown in Figure fTsl 

Angular and projected distances of regions from /3 Pic: (1) 2!'0-3"5 (39-67 AU); (2) 3!'7-5"6 (71-108 AU); 
(3) 6"6-10"0 (127-193 AU); (4) 10;'0-13;'4 (193-258 AU). 

■^Indices in region 4 of the PSF-deconvolved images are suspect because the efficacy of the Lucy-Richardson 
algorithm is diminished in regions of low S/N. 



TABLE 4 

Power-law fits to surface-brightness profiles of component disks" 



Extension 


Component 


Rej 


;ion 




Power-law index (af' 




(arcsec) 


(AU) 


F435W 


F606W 


F814W 


Noilheast 


Main 


4.1-5.6 


80-108 


1.23 ±0.06 


1.42 ±0.07 


1.19±0.08 




Main 


6.6-7.8 


127-150 


5.09±0.12 


4.52±0.11 


4.38±0.25 




Secondary 


4.1-7.8 


80-150 


4.46 ±0.14 


4.63 ±0.14 


4.97 ±0.24 


Southwest 


Main 


4.1-5.6 


80-108 


1.50±0.05 


1.64 ±0.06 


1.58±0.06 




Main 


6.6-7.8 


127-150 


5.23±0.16 


5.18±0.14 


4.96±0.12 




Secondary 


4.1-7.8 


80-150 


3.80±0.10 


3.72±0.13 


3.93±0.10 



"Profiles are extracted from the PSF-deconvolved images and shown in Figure fT^ 
''Seenote(a)ofTable|3l 
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Fig. 1. — 29" x 10" sections of the F606W (Broad V) images of /? Pic (top and middle panels) and the reference star a Pic (bottom panel) obtained with 
the ACS HRC coronagraph. The images are displayed with logarithmic scaling, but without correction of geometric distortion. The dust disk around (3 Pic, 
which is viewed nearly edge-on, is evident without subtraction of the stellar point-spread function. The apparent position angles of the disk in the top and middle 
panels change in accordance with the HSTmW offset of 9°.7. The linear feature seen in all three panels is a component of the coronagraphic PSF whose origin is 
presently unknown. It is coUinear with the HRC's occulting finger, which lies beyond the field of view of each panel. 
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lDg^Q[Surtace Brightness (relative 1o [3 Pic)] 



-10 -9^ -9 -8.5 -B -7.5 -7 ^6.5 




Projected liorizontal distance (AU) 

Fig. 2. — Multiband HRC images of the disk around /3 Pic after subtraction of the stellar PSF, but before deconvolution of the "off-spot" PSF. The panels show 
21'.'2 X 8"3 sections of the F435W (B). F606W (Broad V), and F814W (Broad /) images. The image sections have been rotated so that the northeast extension 
of the disk appears to the left of each panel and the midplane of the outer disk (radius > 100 AU) is horizontal. Each color-coded panel shows the logarithm of 
the disk's surface brightness relative to the star's brightness in that bandpass. The in'egular, blackened regions near the center of each panel reflect imperfect PSF 
subtraction around the occulting spot, especially along the direction of the occulting finger, which is beyond the exhibited FOV. 
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Surface brighlness (relative 1o spine o1 disl<} 
0.2 0.4 0.6 f 





■"2 " l^ 1 T i li I I 

-200 



Projected liorizontal distance (AU) 



Fig. 3. — Same as Figure|2lwith the vertical scal e exp anded by a factor of four. The color-coded images represent the surface brightnesses of the disk relative 
to those measured along the spine of the disk. (See i|2.2l for details.) The expanded vertical scale exaggerates the warp in the inner region of the disk first reported 
bv lBurrows et alj|l995l) . The red dot seen to the left of the "NE" label in the bottom panel is a very red background source located 1 1 !'5 from /3 Pic at a position 
angle of 32.0°. 
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Fig. 4. — Same images shown in Figure l2l after Lucy-Ricliardson deconvolution of the "off-spot" PSF. Pixels lying within a radius of l!'5 (~ 30 AU) of the 
image center (i.e., the location of the occulted star) have been masked and excluded from the deconvolution. Because the Lucy-Richardson algorithm forces 
all pixels to have positive values, the deconvolved images exhibit no negative PSF-subtraction residuals and enhanced, correlated noise at faint signal levels. 
Consequently, photometry of the disk in these regions of the disk is less reliable than in regions having large S/N ratios and small subtraction residuals. 
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Fig. 5. — Same images shown in Figure l3l after Lucy-Ricliardson deconvolution of the "off-spot" PSF. The color-coded isophotes are more similar among 
these deconvolved filter images than among the convolved images shown in Figurel3l Side effects of the deconvolution process are described in FigureBI 
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Fig. 6. — Isophotal map of our F606W images of the disk before PSF deconvolution (top panel) and after PSF deconvolution (bottom panel). Tlie images 
were smoothed with a 9 x 9-pixel boxcar. The interval between isophotes is 1 mag arcsec"^ ; the outermost isophote in each panel marks a surface brightness 
of 20 mag arcsec"^. A circular mask of radius l!'5 (~ 30 AU) has been imposed on the innermost region of the disk to reduce confusion f rom PS F-subtraction 
residuals. Each dashed line represents the least-squares fit to the midplane of the corresponding extension of the inner, secondary disk (see il3.1.1l . 
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Fig. 7. — Ratios of HRC filter images after and before deconvolution of the off-spot coronagrapliic PSF (i.e., images in Figure |4]divided by corresponding 
images in Figure|2j. Tliese images accentuate tlie sharply-peaked midplane of the disk and support the notion that the inner warp is a secondary disk, distinct 
from the main outer disk and inclined from it by ~ 5°. A circular mask of radius 1"5 (~ 30 AU) has been imposed on the innermost region of the disk to reduce 
confusion from PSF-subtraction residuals. 
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Fig. 8. — The vertical scattered-light profiles (i.e., the scattered-light distributions perpendicular to the midplane) of each extension of the disk at a projected 
horizontal distance of 100 AU from /3 Pic. The jagged solid curves are the observed profiles extracted from our F606W images after PSF-deconvolution. Each 
smooth solid curve is the sum of two hybrid-Lorentzian functions (shown as dashed curves; see i|3.1.U that best fits the respective observed profile. The dashed 
curves represent the contributions of the main and secondary disks to the composite vertical profile. The asymmetry of each composite profile and the reversal of 
this asymmetry in opposing extensions of the disk reflect the "butterfly asymmetry" characterized bv .Kalas & Jewitl.il995il . 
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Fig. 9. — Traces of the two "hybrid-Lorentzian" components of the composite vertical profiles (e.g., Figurelsl at horizontal distances of 30-250 AU from /3 Pic. 
The vertical dotted Hnes bound the regions where the composite profiles can be accurately decomposed into two profiles associated with the main (horizontal) 
and secondary (tilted) disks. The solid lines are least-squares fits to the traces along each disk extension in the regions between 80 and 250 AU for the main disk 
and between the vertical dotted lines for the secondary disk. Profile decomposition breaks down within ~ 80 AU of the star, as evidenced by the large deviations 
of the main disk components from their respective linear fits. The noisy traces of the secondary disk within ~ 80 AU have been omitted for clarity's sake. 
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Fig. 10. — Magnified traces of the main disk component observed in our F606W image after PSF deconvolution. The vertical dotted lines mark the 80 AU 
boundary within which the composite vertical scattered-light profile cannot be credibly decomposed into its main and secondary disk components. The solid lines 
are the least-squares fits to the traces along each extension of the main disk at distances of 80-250 AU from {3 Pic. The opposite slopes of the linear fits reflect 
a persistent but diminished "wing-tilt asymmetry," first noted by Kalas & Jewitt ( 1995) at lai'ge distances from the star. The general conformity of the northeast 
extension to the Hnear fit is not observed in the southwest extension. The opposing traces obtained from the F435W and F814W images show similar asymmetry; 
their average vertical positions deviate from the F606W traces by < 0.5 AU (i.e., less than one HRC pixel) beyond 80 AU from the star. 
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Fig. 11. — The projected vertical width of the composite disk as a function of horizontal distance from /3 Pic. The curves trace the full widths of each extension 
of the disk at half {top panel) and one tenth {bottom panel) of the maximum (midplane) brightness. The upper pair of thin soUd and dashed curves in each panel 
show the widths of the northeast and southwest extensions, respectively, measured from the F606W image before deconvolution of the off-spot PSF. The lower 
pair of thin curves show the reduced widths after PSF deconvolution. The short, thick pair of curves in each panel show the widths of the main component of the 
disk between 80 and 150 AU after decomposition and extraction of the tilted secondary disk. All curves have been smoothed with a 9-pixel boxcar. 
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Fig. 12. — Isophotal map of the inner region of the disk obtained from our PSF-deconvolved F606W image. The data have been smoothed with a 3 x 3-pixel 
boxcar. The isophotes represent surface brightnesses of 10-15.4 mag arcsec"^ at intervals of 0.2 mag arcsec"". The dashed lines are the best linear fits to the 
spine of the secondary disk. The locations of the dia metrically-o pposed clumps of 1 8 /xm emission are marked with the letters assigned to them by Wahhai et ^ 
120031) . The bright clump of emission reported bv lTelesco et al J f2005l) is nearly coincident with clump C. No such clumping is evident in scattered optical light. 
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Fig. 13. — Multiband surface brightness profiles measu red a long the spine of each disk extension before PSF deconvolution. The dotted curves are ilcr en'or 
profiles derived from the total-error maps of the images ( i|2.2l and ~ 3% uncertainty in the photometric calibration. The dashed lines are least-squares fits to 
the logarithmic data at angular distances of 2"-3"5, 3!'7-5"6, 6!'-10", and 10"-13!'4 from f) Pic. The slopes of these lines (i.e., the indices of the equivalent 
power-law fits) are given in Tablel^ The vertical lines at ~ 1"5 mark the inner limit of credible photometry. 
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Fig. 14. — Multiband surface brightness profiles measured along the spine of each disk extensio n aft er PSF deconvolution. The dotted curves are ilcr error 
profiles obtained by combining in quadrature the ±lcr errors of the PSF-convolved images (Figure [T3l and local estimates of the corr elated noise imparted by 
the Lucy-Richardson algorithm. The dashed lines are least-squares fits to the logarithmic data for the regions described in Figure [Tsl their slopes are given in 
Tablel3l The vertical lines at ~ 1"5 mark the boundary of the circular pixel mask used during PSF deconvolution. 
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Fig. 15. — Multiband surface brightness profiles measured along the spine of each disk extension after PSF deconvolution. The data have been smoothed with 
a 5-pixel boxcar to improve clarity. The dashed lines at 30 AU mark the boundary of the circular pixel mask used during PSF deconvolution. The uncertainty in 
the photometric calibration of each bandpass is ~ 0.03 mag. 
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Fig. 16. — Muldband surface brightness profiles of the main (top cun'es) and secondary (bottom curves) disk components, obtained from the PSF-deconvolved 
images. The uns moothe d curves trace the maxima of the hybrid-Lorentzian profiles that best fit the vertical profiles of the two components between 80 and 150 AU 
from /3 Pic (see i|3.1.1i . The uncertainty in the photometric calibration of each bandpass is ~ 0.03 mag. 
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Fig. 17. — F606W/F435W and F814W/F435W ratio images of the composite disk after PSF deconvolution. Tlie images have been smoothed with a 7 x 7 
pixel boxcar. The flux ratios are measured relative to those of the star in the respective bandpasses. Within 150 AU of the star, the uncertainties in the flux ratios 
range from ~ 5-10% along the spine of the disk to ~ 25% at projected distances of ±50 AU from the spine. The F814W/F606W image is not shown, but is 
qualitatively similar to the F606W/F435W and F814W/F435W images. 
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Fig. 18. — F435W-F606W and F435 W-F8 1 4W colors of the composite disk relative to those of /3 Pic, obtained before (dotted curves) and after (solid curves) 
PSF deconvolution. The disk's colors were measured along the spines of each extension. The ratios of the filter images were smoothed with a 7 x 7-pixel boxcar 
for improved clarity. The uncertainties of both colors before PSF deconvolution are ~ 3% at 40-150 AU and ~ 8% at 150-250 AU, whereas those after PSF 
deconvolution are ~ 8% at 40-150 AU and ~ 23% at 150-250 AU. 
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Fig. 19. — F435W-F606W and F435W-F814W colors of the northeast (solid curves) and southwest (dashed curves) extensions of the main (thick curves) 
and secondaiy (thin curves) disks, relative to the colors of /3 Pic. The colors are determined from the PSF-deconvolved images using the maxima of the hybrid- 
Lorentzian profiles fitted to the main and secondary disks 80-150 AU from the star ( il3.1.1i . The data have been smoothed with a 3-pixel boxcar for clarity. The 
RMS deviations of the colors in the main and secondary disks are 0.02 mag and 0.1-0.2 mag, respectively. 
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Fig. 20. — (Qeff) versus grain radius for F435W (dotted cun'es), F606W (solid curves), F814W (dashed curves) and five combinations of porosity and 
composit ion: (a) compact grains of pure astronomical silicate ("astrosil," Draine & Lee 1984); (b) compact grains with equal amo unts of astrosil and graphite 
IDraine & Leenl984D: (c) 33% porous grains with equal numbers of astrosil and graphite inclusions tVoshchinnik ov et al.il2005l): (d) 60% porous grains of 
pure astrosil IWolffet alii 998): and (e) 90% porous grains with equal numbers of astrosil and graphite inclusions iVoshchinnikov et al. 2005). The curves are 
computed for an assumed grain number density of n{r) « r^^. 
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Fig. 2 1 . — Simulated ACS F435W-F606W colors (dashed curves) and F435W-F814W colors (solid curves) of the disk (relative to those of /3 Pic) as func tions 
of minimum grain size for various grain compositions and porosities. The colors are derived from the coiresponding values of (2cff) shown in Figure |20| and a 
grain-size distribution of dn « ir^-^ da iDohnanvi 1969). They pertain to the non-icy grains expected within the ice sublimation zone within ~ 100 AU o f the 
star iPantin et al. 1997; Li & Greenberg 1998). The colors of the composite disk in this region, measured from our PSF-deconvolved images (Figure [Tsl . are 
F435W-F606W ^ 0.1 and F435W-F814W !^ 0.2. 



